I 


NASA Contractor Report 3780 


LOAN^OptR 
AF^lOTGI^i^ 
KIRTLAND 



NASA 

CR 

3780 

c.l 


Computer Program for Prediction 
of the Deposition of Material 
Released From Fixed and 
Rotary Wing Aircraft 



Milton E. Teske 


CONTRACT NAS1-16031 
MARCH 1984 




TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


NASA Contractor Report 3780 


□□ 


(=•3338 


Computer Program for Prediction 
of the Deposition of Material 
Released From Fixed and 
Rotary Wing Aircraft 


Milton E. Teske 

Continuum Dynamics , Inc. 
Princeton , New Jersey 


Prepared for 

Langley Research Center 

under Contract NAS 1-1 6031 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Office 


1984 



SUMMARY 


The AGDISP code computes the motion of agricultural materials 
released from aircraft , predicting the mean position of the 
material and the position variance about the mean as a result of 
turbulent fluctuations. This document is the User Manual to the 
Mod 2.0 version of the AGDISP code. 

This manual consists of nine sections. In Section 1 the dis- 
persal dynamics formulations are reviewed, and a summary of the 
available background mean velocity and turbulence fields is given. 
In Section 2 the structure of the program, and of its companion 
plotting code AGPLOT and equivalent profile distribution code 
AGLINE, is highlighted. Section 3 contains a detailed explanation 
of input requirements and code outputs for AGDISP, while Section 4 
contains similar information for AGPLOT, and Section 5, for AGLINE. 
Section 6 examines a test problem: generating the initial condi- 
tions, constructing the input deck, running AGDISP, plotting the 
particle trajectories with AGPLOT, and forming the equivalent dis- 
tribution with AGLINE. In Section 7 a program flow chart is 
discussed, followed in Section 8 by a subroutine-by- subroutine 
description of the code. Section 9 contains details of the sample 
test cases, and the three Appendices identify variables used in the 
code. 
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1. AGDISP FO EMULATION , OPTIONS AND NOMENCLATURE 


A Lagrangian approach is used to develop the equations of 
motion of discrete particles released from aircraft, and a pre- 
dictor-corrector solution scheme is used to solve the resulting 
set of ordinary differential equations. The formulation of the 
particle equations , particularly the ensemble -averaging approach 
used to develop turbulent properties of the particle-atmospheric 
interaction, is detailed in Ref. 1. This section of the User 
Manual examines those variables and nomenclature that are a part 
of the AGDISP code results. 

The aircraft is assumed to be in level flight near the sur- 
face, releasing particles into a (y,z) plane normal to the 
flight direction, x , as shown in Figure 1. The AGDISP code 
solves for the transport of these particles within this plane, 
until they deposit upon the surface or are carried aloft by wind 
or vortex motions. Accuracy of predictions requires accurate 
description of background environment (mean winds and turbulence) . 
The need of accuracy must, however, be balanced by the require- 
ment of computational speed and ease of calculation. Simulations 
need to be performed simply and repeatedly under a variety of 
background conditions , with minimum computer demands and setup 
requirements. When the need is justified, the solution for the 
background environment may be found by a larger, more detailed 
code, such as the NASA WAKE code discussed in Ref. 2. The mean 
velocity and turbulence profiles predicted by these codes may 
then be read by AGDISP to predict particle behavior. As an 
alternate approach the AGDISP code has been fitted with extensive 
simplified flow field options to permit the prediction of part- 
icle motion in a wide range of idealized background conditions. 

The detailed derivations of these options are included in the 
Technical Report (Ref. 1). In this section of the User Manual, 
the nomenclature needed to implement these options is discussed. 

The released particles are assumed to be spherically shaped. 
The particle flight path as a function of time after particle re- 
lease is computed as the locations (y,z) for each particle 
included in the simulation. The velocity of the particle is de- 
noted by (V,W) . Gravity is downward. The particle is affected 
by the background environment; however, the environment is assumed 
to be unaffected by the presence of the particle. The inter- 
action of the particle with the turbulence in the environment 
creates turbulent correlation functions for the particle position 
and the particle velocity, <yv> and <zw> ; for the particle 
velocity variance, <w> and <ww> ; and for the particle posi- 
tion variance, < yy > and <zz> . The square root of these last 
two variables gives the horizontal and vertical standard devia- 
tions of the particle motion about the mean described by y and 
z . 
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In the AGDISP code the particles are assumed to be influenced 
by three forces: 

1. aerodynamic drag; 

2. forces as a consequence of evaporation; and 

3. gravity. 

Inherent in the particle equations is a particle relaxation 
time xp which is essentially the e-folding time for the part- 
icle to come up to speed with the background velocity. The nu- 
merical solution of the particle equations requires the simulation 
time step to be a fraction of Tp . Likewise, the evaporation 
model possesses a time scale Tg during which the particle will 
e-fold in size due to evaporation. The evaporation model (Ref. 3) 
requires that a wet-bulb depression be user-specified. Fluctua- 
tions in evaporation rate are neglected in the AGDISP code. 
Aerodynamic drag is accounted for by using the empirical drag law 
of Langmuir and Blodgett (Ref. 4) . 

The behavior of the particles is intimately connected to the 
local background mean velocity and turbulence fields through which 
the particles are transported. The flow field options available 
in Mod 2.0 are the following: 

Fixed-Wing Fully Rolled Up Tip Vortices 

If roll up of trailing vorticity can be approximated as 
occurring immediately downstream of the wing, the mean velocity 
field that results may be simply characterized by the aircraft 
semispan, circulation and load distribution (either rectangular 
or triangular) . 

Fixed-Wing Betz Roll Up 

Time integration of trailing vorticity rolling up into mult- 
iple discrete vortices and the influence of the increasing vortex 
strength and of the unrolled vortex sheet are available. The 
approach is based on a methodology extended from Betz (Refs. 5 
and 6) and involves specifying the wing load distribution. 

Propeller 

The propeller model is an adaptation of actuator-disk theory 
(Ref. 7) and involves knowing the basic propeller characteristics 
of engine efficiency, shaft rpm, blade radius, aircraft drag co- 
efficient, and aircraft wing planform area. Turbulence effects 
due to the propeller slipstream are included. 

Helicopter in Forward Advance 

The helicopter flow field is modeled partly as a skewed 
hover wake, following actuator-disk theory, and partly as a 
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fully-rolled up rectangular vortex pair. The advance ratio y 
(y > 0 for forward flight) is used to partition the helicopter 
weight between the two flow fields (see Ref. 1 for details). 

Mean Background Crosswind 

Atmospheric crosswind in the (y,z) plane is parameterized 
as a neutral logarithmic velocity profile dependent upon knowing 
the crosswind at a height and the surface roughness (generally 
characterized as 1/30 of the physical height of the surface 
covering) . 

Plant Canopy 

When a plant canopy, such as a stand of fir trees, comprises 
the surface features , the plant areal density of the canopy may 
be entered into the AGDISP code. Modifications to surface rough- 
ness, turbulence and particle paths are made consistent with the 
canopy model of Ref. 8. 

Sup erequi librium Turbulence 

Local background turbulence level may be simply specified by 
specifying the crosswind velocity, by fixing a constant value 
(typically one percent of the mean wind speed) , or by optioning 
the AGDISP code to compute the turbulence field consistent with 
the given mean velocity field. This option solves the nonlinear, 
algebraic equations of superequilibrium turbulence (Ref. 9) , and 
typically requires a doubling of the computer time needed for 
solution of the particle means and standard deviations. 

Terrain 


A nonhorizontal, level surface may be assumed. 
WAKE Results 


The AGDISP code has been fitted with the option of entering 
any externally-computed mean velocity and turbulence field, as 
long as the computer-generated plot file corresponds to the for- 
mat specified in Section 3 of this manual. 
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2. AGDISP ARCHITECTURE 


The AGDISP code consists of three programs: AGDISP for estab- 
lishing the desired background fields and computing the particle 
trajectories, AGPLOT for plotting the resulting particle solutions 
and AGLINE for constructing an equivalent Gaussian profile distri- 
bution from the multiple-particle solution generated by AGDISP. 

These programs are configured as interactive programs reading 
their needed input data from disk files and writing output to the 
terminal. AGDISP reads an input file of data cards (valid cards 
are discussed in Section 3) , types information to the terminal as 
the run proceeds, generates printed output in a separate file, 
stores particle trajectory information in a plot file for subse- 
quent plotting, and reads WAKE plot file results from a separate 
data file (if invoked). AGPLOT reads the desired particle plot 
file and interactively questions the user to supply necessary scal- 
ing data for the plot construction. AGLINE reads the particle plot 
file to construct an equivalent profile distribution (compatible 
for entry into the FSCBG Forest Service code, Ref. 10), generating 
printed output and a Gaussian pattern plot. These programs have 
straight-line flow charts (examined in Section 7 of this manual; 
an overview of their construction follows in this section of the 
User Manual, including discussion of accuracy, program limitations, 
and warning messages. 

The AGDISP code requires the data entry explained in Section 

3. This input data is read and processed, causing termination of 
the program if errors are discovered or needed cards are missing. 
Initialization follows: particle locations and velocities, vortex 
strength, Betz roll up, plant canopy, WAKE plot file. The part- 
icle equations are then integrated repeatedly until one of two 
termination conditions occur : 

1. maximum simulation time is reached; or 

2. all of the particles deposit on the surface. 

Each integration step in turn includes the following: 

1. updating the background mean velocity and turbulence 
fields by the WAKE plot file, Betz roll up, vortex 
motion, plant canopy; 

2. equation system predictor of the particle equations; 

3. two equation system correctors of the particle equations 
and 

4. incremental file save of the particle results for sub- 
sequent plotting by AGPLOT and use by AGLINE. 

Solution of the particle equations is by a modified predictor- 
corrector scheme (Ref. 11) employing two correctors to enhance 
solution accuracy. A fixed time step size is determined by the 
program knowing the particle and evaporation time scales ip and 
TE . The algorithm involved was selected to give accurate results 
to the particle equations without impacting computer time involved. 
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Nonetheless, AGD1SP is flexible enough to admit a wide 
variety of input conditions, some of which may demand increased 
computer time. For these cases, "$$$ WARNING” messages will 
flag such events, to warn the user of cases where input conditions 
may require more computer time. The messages that may appear 
are the following: 

1. SMALL PARTICLE INVOKED occurs whenever the simulation 
follows particles of 50 microns or less. In this case, 
the particle time scale can become small, and the com- 
puting needs, substantial. 

2. SUPEREQUILIBRIUM TURBULENCE INVOKED occurs in those 
instances where the desired turbulence background 
solves the superequilibrium equations. This option 
appears to be most severe, in general doubling the 
computer needs of the AGDISP code. 

3. EVAPORATION INVOKED occurs because the equation time 
step must be adjusted downward as the particles become 
smaller to account for evaporation effects. The com- 
ments of message 1 then apply. 

4. MANY PARTICLES INVOKED occurs whenever the number of 
particles to be followed exceeds ten. Ten is an arbi- 
trarily chosen number, perhaps too low, but serves to 
warn the user that program solution is linear; i.e., 
twenty particles take twice as long as ten particles. 

5. LONG SIMULATION TIME INVOKED occurs whenever the maxi- 
mum expected simulation time normalized by particle 
size exceeds a given level. This level may also be 
low, but it nonetheless serves the same purpose as the 
level in message 4. 

6. WAKE PLOT FILE INVOKED is written when this option is 
used because the overhead associated with reading an 
externally constructed data file on the Langley computer 
system may vary greatly with the size of the data file. 

The limits of the AGDISP code are the following: 

1. All of the idealized background mean velocity and turbu- 
lence fields are for neutral environments. 

2. The total number of particles that can be included in 
the simulation cannot exceed 60. The error message 
INCORRECT NUMBER OF PARTICLES will be invoked and AGDISP 
will terminate. 

3. The Betz roll up cannot enter more than 100 discrete 
circulation values along the wing (T as a function of 
y). The error message ERROR IN CIRCULATION DATA INPUT 
will be invoked and AGDISP will terminate. 
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4. The Betz roll up procedure cannot handle more than four 
discrete vortices rolling up on each wing. In this 
case, the error message BETZ WILL ROLL UP MORE THAN 4 
VORTICES will be invoked and AGDISP will terminate. 

5. The plant canopy input data cannot permit more than 100 
discrete entries to define the canopy vertical profile 
shape. The error message ERROR IN PLANT AREA DENSITY 
INPUT will be invoked and AGDISP will terminate. 

The AGPLOT code interactively plots the AGDISP plot file 
contents on Tektronix 40xx terminals using the basic CDC plotting 
commands. Input is the desired plot file; after scanning this 
file AGPLOT questions the user about particle plot character- 
istics and sizes of axes, then plots the selected data. Following 
this plot, the ground deposition pattern may also be plotted, if 
any particles have hit the surface. The length of time needed to 
plot the results depends principally upon the number of particles 
being plotted and the amount of data to plot (which in turn 
depends on some of the inputs to AGDISP). 

AGLINE quizzes the AGDISP plot file to compute the equivalent 
Gaussian distribution of the multiple-particle solution. Since 
the position and turbulent standard deviation of every particle 
in the simulation is computed by AGDISP as a function of time after 
particle release, the equivalent mean position and standard devi- 
ation may be determined by appropriate summation and integration 
at each time over all particles in the cross-section. A measure 
of the compatibility of the Gaussian distribution replacing the 
multiple-particle Gaussian distributions is made in the AGLINE 
code by computing a figure of merit ranging from 0 to 1. When 
the equivalent Gaussian no where represents the multiple-particle 
distributions, the figure of merit equals zero. When the equiva- 
lent Gaussian is every where identical to the multiple distribu- 
tion the figure of merit is unity. The user then ultimately 
judges the compatibility of fit and either uses the "best fit" 

(the highest figure of merit) or reruns AGDISP with additional 
detail to improve accuracy. To aid in this compatibility assess- 
ment, a plotting package similar to AGPLOT is included in the 
AGLINE code to plot the multiple-particle distributions and the 
equivalent Gaussian pattern, one upon the other. 
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3. AGDISP INPUTS 


This section of the AGDISP User Manual details the input 
cards to the code. All data entry is in free format, wherein 
card data is separated by commas or blank spaces. This conveni- 
ence offers ease of formatting the data but requires that every 
data card used have all of its data elements present, even if 
they are zero. Unless noted below as integer values, data is 
entered as real numbers (with decimal points). The MKS system 
is used throughout. 

A typical interactive procedure file for the AGDISP code is 
illustrated in Table 1. Here data must be present in a reader 
file attached to unit 4, interactive output is written to the 
terminal on unit 6, printed output is written to a file on unit 9, 
particle trajectory data results are written to a data output file 
on unit 8 (for later plotting by AGPLOT) , and WAKE plot file data 
is entered on unit 10. All card input data is echoed on the ter- 
minal, improper data is flagged, and missing data cards are ex- 
posed. All data cards begin with a four-digit identification in 
columns 1-4, with the rest of the data in free format. The order 
of the cards is important and must follow ever-increasing identi- 
fication numbers. The AGDISP code has been programmed to verify 
this order. In addition, certain available options are inconsis- 
tent with each other. The AGDISP code has been programmed to 
trap these inconsistencies. Unsupported cards are flagged with 
the message AGDISP CODE DOES NOT SUPPORT CARD NUMBER, while mis- 
sing data cards are flagged with INPUT DOES NOT FULLY INITIALIZE 
AGDISP RUN. With free format there exists the chance that all 
necessary data does not appear on the appropriate card; in this 
case, the AGDISP code will run out of data cards before it ration- 
alizes all of its data pointers. The error message will read 
INSUFFICIENT DATA BEFORE CARD following the card where data ended. 
This error message is perhaps the most devious because it may 
require a detailed examination of all data cards to determine the 
missing data point. 

On the other hand, because all of the data is echoed on the 
terminal, error messages should be traceable in a systematic 
manner. A careful check should be made the first time a new case 
is started, and at least until the user is confident of the pro- 
gram input requirements . Only by looking at what the program 
thinks is inputted will the user be able to verify that what was 
inputted was correct. In almost all cases, the code makes no 
check of input validity, either signs or magnitudes. Gravity 
is built into the program as 9.8 m/sec^. Every AGDISP run re- 
quires the entry, at least, of cards 0010, 0020, 0050 and 0060. 

Detailed description of the available data cards follows, 
including the special error messages they invoke. Examples of 
all cards are shown in Table 2. 
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Card 0000 


is a comment card, any number of which may be placed 
anywhere in the input deck. The 76 columns of com- 
ment are merely reproduced on the terminal. 

Card 0010 is the program card requiring two entries. The first 
is the maximum time in seconds for the simulation to 
run, a time which will be shortened if all of the 
particles impact the surface before this time is 
reached. The other entry is the simulation plane 
entry (integer) configuring either a half -plane 
(y > 0) solution (a value of 1) or a full plane 
(_oo < y < +co) solution (a value of 2) . The full plane 
solution should be invoked if a crosswind, propeller or 
terrain exists. Single particle release without an air- 
craft should use the full plane value of 2. Card 0010 
is a required input card. 


Card 0015 is the terrain slope angle card requiring the entry 

of local ground slope in degrees. A positive terrain 
angle raises the right half of the surface as viewed 
in Figure 1, with the origin remaining along the 
initial aircraft vertical centerline. The terrain is 
characterized by a locally straight surface so that 
all of the simplified flow field options in AGDISP 
remain available. A full plane solution is necessary. 
Canopy and crosswind effects remain parallel to the 
tilted surface. This feature cannot be invoked with 
the WAKE plot file option. Card 0015 is an optional 
card. 

Card 0020 is the aircraft characteristics card requiring six 

entries. The first is the mean velocity flag (inte- 
ger), taking one of six values: 

3 denotes helicopter entry; 

2 denotes a rectangularly loaded, fully 

rolled up tip vortex; 

1 denotes a triangularly loaded, fully 
rolled up tip vortex; 

0 denotes Betz roll up from a given 
circulation pattern; 

-1 denotes WAKE plot file entry 
(explained under card 0050); 

-2 denotes a nonaircraft run, appropriate 
for single particle release. 

An entry of 3 requires a 0030 card; entry of 2 or 1 
requires a 0022 card; and entry of 0 requires multiple 
entries of 0025 cards. The second (integer) input is 
the crosswind flag. A value of 0 implies no crosswind; 
a value of 1 does, and requires a 0028 card. The third 
entry is the semispan of the aircraft in meters (this 
is also the rotor radius for a helicopter, the initial 
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y vortex coordinate for flag = 2 and twice the initial 
y coordinate for flag =1). The fourth entry is the 
nominal height of the aircraft wing above the surface 
in meters (this distance becomes the release height 
for single particle release, and is the assumed 
nominal release height for multiple particles). For 
fully rolled up tip vortices , this height is the 
initial z coordinate of the vortex centerline. For 
Betz roll up data this height is the z coordinate 
of the initial vortex sheet (assumed horizontal). 


The fifth entry is the flight speed of the aircraft 
in m/sec. The sixth entry is the biplane configura- 
tion flag (integer) and is 0 for a single wing air- 
plane and 1 for a biplane (with further information 
supplied on a 0021 card). Card 0020 is a required 
input card. 


Card 0021 requires entry on card 0020, and specifies the bi- 
plane characteristics with three entries. The first 
is the vertical distance in meters from the main wing 
location specified on card 0020 to the biplane wing. 
The second entry is the semispan of the biplane wing 
entered as a fraction of the semispan of the main 
wing (if the wings were equal in length this entry 
would be 1.0). The third entry is the vortex strength 
of the biplane wing entered as a fraction of the vor- 
tex strength of the main wing (for equal strengths 
this entry would be 1.0). 


Card 0022 requires entry on card 0020, and specifies the rect- 
angular or triangular loading circulation value, 
m^/sec, for a fully rolled up tip vortex. 


Card 0025 requires entry on card 0020, and enters the Betz wing 
load distribution; each card consists of two entries. 

The first entry is the position in meters measured 
from the wing root monotonically towards the tip, 
while the second entry gives the corresponding circu- 
lation value in m^/sec at this location. The final 
card is signaled by entering the last position as the 
negative of its actual value. This card then forces 
initialization of the Betz roll up procedure. This 
initialization (and the ensuing roll up) invokes extra 
printer output summarizing the roll up process. One 
extra, rather specific error message will be invoked 
if the input load distribution is not "smooth" in a 
sense determined by the code. In this case the error 
NONDISCRETE DISTRIBUTION LOCATION is displayed, sug- 
gesting that the user inspect the generated circula- 
tion derivatives for lack of monotonicity. 

As the Betz roll up continues , the particles being 
tracked will be influenced by every- increasing strength 
vortices , whose positions will approach those of fully- 
rolled up vortices. The AGDISP code includes the effect 
of the unrolled up sheets on the motion of the particles. 
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Card 0028 requires entry on card 0020, and is the crosswind 

card entering the three data values needed to esta- 
blish the neutral crosswind velocity profile shape. 
The crosswind velocity is the horizontal velocity in 
the y direction. The first entry is the mean wind 
velocity in m/sec at the altitude given by the second 
entry, in meters, while the third entry gives the sur 
face roughness height z Q in meters. 

Card 0030 requires entry on card 0020, and is the helicopter 
card consisting of two entries. The first is the 
weight of the helicopter in N , and the second is 
the forward advance ratio y , a value greater than 
zero. The helicopter flow field is idealized as a 
propeller downwash field and a rectangularly loaded, 
fully rolled up vortex pair. 

Card 0040 is the propeller data card consisting of six entries. 

The first is the drag coefficient of the airplane, 
while the second is the planform area of the airplane 
in m2. The third entry is the aircraft propeller 
efficiency and the fourth is the shaft rpm. The 
fifth is the propeller blade radius in meters and the 
sixth entry is the incremental distance in meters of 
the shaft centerline above or below the nominal 
release height given on card 0020 (the propeller is 
assumed to be at the airplane centerline, y = 0) . 

Card 0050 is the turbulence data card containing three entries. 

The first is an integer flag carrying one of the 
following values: 

-1 invokes superequilibrium turbulence; 

0 assumes a fixed value of turbulence given 
on this card; 

1 \ 

r , specifies the turbulent components in the 
’ attached WAKE plot file (invoked with a 

3' -1 entry on card 0020). 


The other two entries on this card are the maximum 
value of the background turbulence q2 in m^/sec2 
(this value is incremented by the presence of a 
crosswind) , and the maximum value of the background 
turbulent macroscale A in meters. The macroscale 


is assumed to be zero at the surface and to grow 
linearly with height up to the maximum value speci- 
fied. If the inputted maximum value of A is suffi- 


ciently large, the macroscale will be proportional 
to z ' (the proportionality factor is 0.65) for all 
particle height positions. Card 0050 is a required 
input card. 
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Depending on the value of the integer flag, the 
AGDISP code may not use either of the other entries 
on this card (they must still be present, of course, 
even if they are set to 0.0). The integer entry on 
card 0050 defines the turbulent background character- 
istics and how AGDISP computes them. The determina- 
tion of the background turbulence should not be 
confused with the use of a WAKE plot file (a -1 entry 
on card 0020). The presence of a WAKE plot file 
implies that the background horizontal and vertical 
velocity components are always found in the WAKE plot 
file. The integer entry on card 0050 defines whether 
additional background components are also present on 
the WAKE plot file. If a WAKE plot file is not in- 
voked (the entry on card 0020 Is something other than 
-1), only the -1 and 0 entries on card 0050 are 
admissible. The -1 entry on this card invokes the 
super equilibrium option, overriding the maximum value 
of q2 entered on this card but using the maximum 
value of A . A0 entry on card 0050 uses both maxi- 
mum values of q2 and A . When a WAKE plot file is 
present, the -1 entry again invokes the superequili- 
brium option and overrides the maximum value of q^ ; 
the 0 entry uses both maximum values (in these in- 
stances, of course, the horizontal and vertical back- 


ground velocities are recovered from the WAKE plot 
file). Entries of 1 , 2 or 3 are valid only when a 
WAKE plot file exists. The 1 entry assumes that the 
WAKE plot file contains not only the horizontal and 


vertical background 
ground turbulence 
value entry for q2 


velocities but also the back- 
q2 . In this case, the maximum 
on card 0050 is overridden but 


the macroscale data is used. For entries of 2 or 3 


the macroscale is expected to be on the WAKE plot 
file, and neither of the maximum values on card 0050 


is used. 


The invoking of the superequilibrium option forces the 
spatial solution of the super equilibrium equations 
using an iterative solution of a linear equation sol- 
ver to determine the algebraic solution to a set of 
nonlinear equations. Two specific errors may arise: 
LINEAR SOLVER ERROR IN SUPEREQUILIBRIUM and ROOT FINDER 
ERROR IN SUPEREQUILIBRIUM. Either error is serious, 
and any code results obtained after generation of the 
error are suspect. 

A -1 entry on card 0020 invokes the WAKE plot file 
and forces extra output to the printer. The WAKE 
plot file must be constructed as a sequential binary 
data file (an example is given in Section 8). Con- 
sistency of WAKE plot file contents and AGDISP input 
data is checked , with errors invoking WAKE PLOT FILE 
VARIABLE ERROR. The WAKE plot file configuration 
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must be one of the following: 

1. If the turbulence is super equilibrium (the 
turbulence flag is -1) , the plot file vari- 
able list must consist of the cross-plane 
velocities V and W ( two variables) ; the 
turbulence components are determined from 
super equilibrium, and the macroscale is com- 
puted algebraically using the data on card 
0050. 

2. If the turbulence is fixed at a specified 
value (the turbulence flag is 0) , the turbu- 
lent components are assumed isotropic, the 
macroscale is computed algebraically, and the 
plot file variable list must contain only V 
and W (two variables) . 


3. If the turbulence flag is 1, the WAKE plot 
file must additionaly contain q^ ( three 
variables), with the turbulent components 
assumed isotropic and the macroscale computed 
algebraically . 

4. If the turbulence flag is 2, the WAKE plot 
file must additionally contain A ( four 
variables) , with the turbulent components 
assumed isotropic. 

5. If the turbulence flag is 3, the WAKE plot 

file contains all of the data needed to 
specify the mean velocity and turbulent back- 
ground profiles, ( six variables): V , W , 

q^ , A , <w> and <ww> . 

The first nineteen words of the WAKE plot file are a 
comment header field. The next integer entry gives 
the number of variables stored in the WAKE plot file 
(2, 3, 4 or 6 as explained above). This number must 
be consistent with what AGDISP expects. Following 
this integer is a header field of length equal to the 
number of variables. The next integer on the WAKE 
plot file is the number of y (or horizontal) mesh 
points in the plot file, followed by these mesh 
values in meters. The next integer is the number of 
z (or vertical) mesh points, followed by these mesh 
values, again in meters. The code restricts these 
entries to 16 or less. Larger values invoke WAKE 
PLOT FILE MESH OUT OF RANGE. 

The profile data follows. Each time slot repeats the 
same pattern. It begins by a single word which is 
the time saved in seconds. A negative time denotes 



end-of data. The data for each variable on the file 
(specified in the order above) follows, by giving all 
of the y values of the first variable at the first 
z position, then all of the y values of the first 
variable at the second z position, on to the last 
z position; then on to all of the y values of the 
second variable at the first z position, etc. , 
until the values of all of the variables at all of 
the y and z positions have been given. The next 
time value follows, repeated to the end of the plot 
file. Velocities are in m/sec, scales in meters, 
and turbulence in m^/sec^. The warning message 
WAKE PLOT FILE EXTRAPOLATION is output the first 
time spatial extrapolation must be used by the AGDXSP 
code during interpolation for the variables contained 
on the WAKE plot file. For times beyond the entries 
in the WAKE plot file, the spatial profiles nearest 
in time to the time being solved in AGDISP will be 
used. 

Before integrating the particle equations , AGDISP 
reads the entire WAKE plot file through to the nega- 
tive end-of -data . If the file is constructed in- 
correctly, the error message PREMATURE END OF WAKE 
PLOT FILE REACHED is invoked, and AGDISP terminates. 

Card 0055 enters the canopy plant area profile; each card con- 
sists of two entries. The first entry is the z 
position in meters beginning at the surface and in- 
creasing monotonically to the top of the canopy. 

The second entry is the value of the plant area 
density (m^/m^) corresponding to the z location. 
The final card 0055 is signaled by entering the last 
z position as the negative of its actual value. 

This card entry then initializes the canopy calcula- 
tion, computing the displacement thickness of the 
canopy, forcing modification to the crosswind velocity 
and turbulence within the canopy. During the AGDISP 
run, the trajectories of the aircraft vortices will be 
altered upon entering the canopy. 

Card 0060 is the particle data card, containing six entries. 

The first entry is the integer number of particles 
to be followed in the half-plane (y > 0); the second 
entry (integer) is equal to 1 if a particle is to be 
released at the centerline of the aircraft (y = 0); 
otherwise 0. Thus, a half -plane solution will track 
total particles equal to the sum of these two inte- 
gers; a full-plane solution will track total particles 
equal to the sum of twice the first integer plus the 
second integer. Additionally, if these two integers 
are entered positive, AGDISP will position the part- 
icles uniformly along the wing (for a first integer 
entry of 1, the particle will be positioned at 1/2 
the semispan; for a first integer entry of 2, the 
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Cards 0061 


Card 0062 


Card 0065 


particles will be positioned at 1/3 and 2/3 the semi- 
span, etc,). If these two integers are entered 
negative, AGDISP will expect sufficient 0061 cards in 
the input deck to initialize all of the particles . 

For both integers equal to zero, AGDISP assumes a 
single particle release at the centerline position. 

The third entry is the vertical position in meters 
off-setting the particle release point from the height 
of the wing given on card 0020. The fourth entry is 
the diameter in microns of the released particles 
(only one particle size can be released during any one 
run of AGDISP) . The fifth entry is the specific 
gravity of the released particle. The sixth entry 
is the (integer) evaporation flag, set to 0 for no 
evaporation or 1 when evaporation is invoked. Evapor- 
ation requires the entry of card 0065. Card 0060 is 
a required input card. 

require entry on card 0060, and are the discrete part- 
icle location entry cards. Each card contains two 
entries. There must be enough 0061 cards to satisfy 
the requirements of the 0060 card. If a full-plane 
solution is invoked all particle locations must be 
specified. If a particle is released along the center 
line, that particle must be the last card 0061 for 
purposes of computing the ground deposition correctly. 
The first entry is the y value in meters along the 
wing, while the second entry is the vertical position 
in meters off-setting the particle release point from 
the height of the wing given on card 0020. The verti- 
cal position off-set on card 0060 is disregarded. 
Discrete particle locations allow for fine tuning of 
nozzle locations , and are necessary if nozzles are not 
uniformly distributed along the wing or if the spray 
boom is not parallel to the surface. 

is the particle initial condition data card, and 
consists of four entries: the initial horizontal 
velocity in m/s ec; the initial vertical velocity in 
m/sec; the initial spatial variance of the particle 
path in m^; and the initial velocity variance of the 
particle in m^/sec^. This card is optional; all 
initial conditions are set to zero unless modified by 
card 0062. 

requires entry on card 0060 , and is the evaporation 
data card, consisting of two entries; the wet bulb 
temperature difference in deg. C, and the particle 
diameter in microns, below which the user considers 
evaporation to have ceased. This latter number is 
important because the smaller the particle becomes , 
the smaller the step size must be to maintain numeri- 
cal accuracy in solving the particle equations. 
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This completes the input discussion to the AGDISP code. The 
calculation is then performed, with periodic output to the printer 
of the particle positions , variances , and turbulent correlations. 
When a particle intersects the surface, it is no longer computed, 
and, therefore, its location is no longer printed. At the end of 
the run, the surface impact statistics are summarized and the 
deposition fraction (the mass fraction of particles released that 
have reached the ground in the desired simulation time) is 
computed. 
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4 . AGPLOT INPUTS 


The AGPLOT program may be invoked to plot the resulting 
trajectories and deposition pattern. The code asks questions of 
the user in the interactive environment, and with this informa- 
tion processes the data and plots the results. Again, entry is 
in free format. A typical procedure file for AGPLOT is shown in 
Table 3. Here the particle trajectory data file (constructed by 
an AGDISP run) is attached to unit 8, interactive output is 
written to the terminal on unit 6 , and a system graphics library 
is accessed to resolve the external plotting subroutine calls. 

All graphics output is displayed on the terminal. The line-by- 
line questions asked by AGPLOT are the following: 

ENTER TERMINAL DESIGNATION requires the entry of -1, 0, 1 or 
2. Entries of -1 and 0 pertain to a Tektronix 4025, with -1 
producing a horizontal plot (the surface z - 0 is horizontal 
on the screen) and 0 producing a vertical plot (the surface z = 0 
is vertical on the right edge of the screen) . The entry of 1 or 2 
pertains to a Tektronix 40 lx and invokes a horizontal plot using 
CDC graphic software (a response of 2 does not produce a hard 
copy). In all cases the plot fills the screen. 

The desired AGDISP data file is then read, and the input 
cards used to generate the results are typed on the screen to 
serve as verification before plotting. 

FULL PLANE PLOT is asked whenever the AGDISP run was a half- 
plane run. Response is YES or NO. No such request is necessary 
when a full-plane calculation is being plotted. 

ENTER Y SCALE MIN, Y SCALE MAX, Y SCALE DELTA asks for the 
horizontal scale information (in meters); at what y position 
the scale starts, at what y position the scale ends, and the 
scale increment at which to place tick marks. The AGPLOT code 
has a scale checking routine and will invoke a warning message 
whenever the scale delta is not an integer fraction of the over- 
all scale size. Also, three errors will be trapped: MAX SCALE 
LESS THAN MIN SCALE: SCALE DELTA LESS THAN ZERO; and TOO MANY 
SCALE DIVISIONS whenever more than ten tick marks are required. 

PLOT PARTICLE MEANS requests YES if the mean position of 
the particles are to be plotted; NO if they are not. The part- 
icle mean trajectories are plotted as solid lines. 

PLOT PARTICLE STANDARD DEVIATIONS requests YES if the tur- 
bulence standard deviations of the particle positions are to be 
plotted about the mean; NO if they are not. If neither means nor 
deviations is requested, the code will move down to the questions 
on ground dispersal pattern. The particle standard deviations 
are plotted as dashed lines. 
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TAG LOCATION AT TIME INCREMENTS requests whether the 
particle mean positions should be tagged (by a +) at uniform time 
increments during the simulation. The size of the + is scaled 
by the size of the particle standard deviation at the (y,z) 
position of the +. If the response is YES, the code will 
further ask ENTER TAG TIME INCREMENT as incremental seconds be- 
tween tags. If the original response is NO, no further informa- 
tion is requested. 

ENTER Z SCALE MIN, Z SCALE MAX, Z SCALE DELTA requests the 
vertical scale data of minimum value, maximum value, and tick 
increment in meters. With the input of this data, all necessary 
information is available to plot the trajectory data. Before 
plotting, though, one further question is asked: 

PLOT SCALE INCREMENTS EQUAL asks whether the y and z 
scales should be plotted in a one-to-one ratio on the terminal. 
This option makes sense only if visual advantage exists in pre- 
senting a more physically correct picture of the trajectory 
pattern. After response to this question, the AGPLOT code will 
plot the desired trajectories and invoke a copy of the screen. 

It will then move on to plotting the ground deposition pattern. 

PLOT GROUND DEPOSITION PATTERN asks whether the deposition 
pattern is desired. If the answer is NO, the AGPLOT code will 
terminate. If YES, the deposition pattern will be computed 
using Gaussians whose magnitudes are consistent with the size of 
the particle on impact, whose center location corresponds to the 
y position of the mean as the particle impacted the surface, and 
whose standard deviation is equal to the standard deviation of 
the particle just before impact. If none of the particles have 
hit the ground, the message NO GROUND DEPOSITION is typed, and 
the AGPLOT code will terminate. Otherwise the code will type: 
ENTER G SCALE MIN, G SCALE MAX, G SCALE DELTA to request the 
deposition scale values of minimum, maximum and increment. After 
this the deposition pattern is plotted and copied. The AGPLOT 
code then terminates . 
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5 . AGLINE INPUTS 


The AGLINE program is invoked by using a procedure file 
similar to the one illustrated in Table 4. Program input data 
comes from terminal response by the user and trajectory data from 
the AGDISP plot file on unit 8. A printer output file is con- 
structed on unit 9 to document the analysis of the plot file and 
all final results. The desired configuration position is plotted 
on the terminal. The AGLINE code requests and actions are the 
following: 

ENTER TERMINAL DESIGNATION requires the entry of -1, 0, 1 
or 2. Entries of -1 and 0 pertain to a Tektronix 4025 (horizontal 
and vertical plots, respectively), while the 1 and 2 entries per- 
tain to Tektronix 401X terminals (a response of 2 does not produce 
a hard copy) . 

The inputted AGDISP plot file is then processed. The input 
data deck used to generate the results are typed on the screen to 
serve as verification before plotting. Equivalent Gaussian cal- 
culations are made for every time on the plot file; the mean posi- 
tion, standard deviation and figure of merit are written to the 
printer output file. All times on the plot file are given a mono- 
tonically increasing sequence number. The critical sequence 
numbers , the time associated with each of them and their respective 
figure of merit are typed to the screen whenever: 

1. the first sequence number is reached; 

2. the last sequence number is reached; 

3. the figure of merit reaches a local minimum; 

4. the figure of merit reaches a local maximum or 

5. a particle first comes within a standard deviation 
of the ground. 

After processing the entire plot file and displaying the above 
summary data, the AGLINE code requests: 

ENTER DESIRED SEQUENCE NUMBER to select the sequence position 
in the plot file at which to display the equivalent Gaussian. If 
a specific figure of merit is desired the appropriate response is 
0 in which case the code will request: 

ENTER DESIRED FIGURE OF MERIT will cause the code to seek 
the first occurrance of the entered figure of merit. If the ent- 
ered value cannot be found the error message DESIRED FIGURE OF 
MERIT NOT FOUND will be invoked and AGLINE will terminate. Other- 
wise the search will be successful and the equivalent Gaussian 
data will be written to the printer output file. If the equivalent 
Gaussian vertical standard deviation places the vertical mean into 
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contact with the surface the warning message GROUND ENCOUNTER BY 
EQUIVALENT GAUSSIAN will be typed to the screen. 

The particle and equivalent Gaussian distributions are now 
plotted. Consistent with the AGPLOT procedures AGLINE will ask: 

ENTER Y SCALE MIN, Y SCALE MAX, Y SCALE DELTA for horizontal 
scale information (in meters) ; at what y position the scale 
starts, at what y position the scale ends, and the scale incre- 
ment at which to place tick marks . 

ENTER Z SCALE MIN, Z SCALE MAX, Z SCALE DELTA for vertical 
scale data of minimum value, maximum value, and tick increments 
in meters. With the input of this data, all necessary informa- 
tion is available to plot the distribution data. Before plot- 
ting, though, one further question is asked: 

PLOT SCALE INCREMENTS EQUAL asks whether the y and z 
scales should be plotted in a one-to-one ratio on the terminal. 
After response to this question, the AGLINE code will plot the 
Gaussian distribution and invoke a copy of the screen. 

The resulting plot represents each multiple-particle distri- 
bution as dashed lines and the single equivalent distribution as 
solid lines. Around each center position normalized contour lines 
are plotted at one- eighth, one- fourth, one-half, one and two times 
the standard deviation. 
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6. EXAMPLE TEST CASE 


This section of the User Manual illustrates the use of the 
AGDISP code by presenting an example test case. The following 
problem is posed (see Figure 2 for a schematic). 

A single-wing aircraft weighing 17300 N and flying at 
35 m/sec is dispersing 150 micron water droplets from four loca- 
tions on a spray boom positioned beneath the wing 3 meters above 
the surface. The wing has a chord of 1.5 meters and a semispan 
of 6 meters, with a slight dihedral such that the rectangular 
wing loading rolls up a tip vortex with a vertical centerline 
0.75 meters above the level of the spray boom. The spray nozzles 
are located at -3.0, 0.5, 2.6 and 4.8 meters measured from the 
airplane centerline. 

Complicating the wake of the airplane is the presence of its 
propeller, located one meter above the spray boom level, with a 
blade radius of 1.5 meters and an efficiency of 0.8 at a rated 
speed of 500 rpms . For convenience, the drag coefficient of the 
airplane will be taken to be 0.1. The airplane is in straight 
and level flight with a slight crosswind blowing from left to 
right at a velocity of 2 m/sec at 5 meters (nearby instrumented 
tower data) . The ground cover will be assumed to have a rough- 
ness of 2 centimeters. 

The above problem definition is sufficient to generate 
initial conditions to the AGDISP code. The input data cards are 
constructed using the documentation in Section 3, and become the 
following. 

The simulation time will be taken to be 10 seconds, enough 
time to allow ground impact of the three inboard particles (the 
fourth particle will probably be "trapped" in the tip vortex and 
carried left-to-right) . The presence of the crosswind and pro- 
peller require a full-plane solution. Thus, card 0010 will read: 

0010 10.0 2 

Rectangular wing loading will be assumed; hence the velocity 
flag is 2, requiring the entry of a 0022 card. A crosswind re- 
quires a 0028 card. The semispan of the airplane is 6 meters; 
the height of the fully rolled up tip vortex is 3.75 meters above 
the surface; the flight speed of the airplane is 35 m/sec; and 
the airplane is a single wing. Thus, card 0020 will read: 

0020 2 1 6.0 3.75 35.0 0 

The rectangularly loaded wing will have a tip circulation 
determined by the formula F = L/2spU where L is the weight of 
the airplane, s is the semispan, p is the air density and U 
is the flight speed. For standard atmospheric air p = 1.2266 kg/m'* 
and T becomes 33.6 m^/sec. Thus, card 0022 will read: 

0022 33.6 
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The crosswind velocity is 2 m/sec at a height of 5 meters, 
with a surface roughness of 0.02 meters. Thus, card 0028 will 
read : 

0028 2.0 5.0 0.02 


The drag coefficient is assumed to be 0.1; the propeller 
efficiency is 0.8 at 500 rpms; and the propeller radius is 
1 . 5 meters . The planform area of the wing is equal to the span 
times the chord length, or 18 m^. The propeller centerline is 
one meter above the spray boom but only 0.25 meters above the 
height given on card 0020. Thus card 0040 will read: 

0040 0.1 18.0 0.8 500.0 1.5 0.25 


The assumption will be made that the crosswind determines 
the background turbulence level, constant for all points in 
space; with the entry of 0.0 for the fixed turbulence level, the 
code will add the level generated by the crosswind. The assump- 
tion will also be made that the turbulent macroscale is linear 
with height; entering a large maximum value will ensure this 
linearity. Thus, card 0050 will read: 

0050 0 0.0 30.0 

Four particles are to be user-specified and released into a 
full-plane simulation. Thus, the first two entries on this card 
are -2 and 0 (no particle is at the centerline). The code will 
then expect four 0061 cards. The vertical release position is 
disregarded because the points are user-specified (if they were 
uniformly placed by the code, the entry would be -0.75 meters to 
recover the 3 meter height of the spray boom) . The particles 
are 150 microns in diameter with a specific gravity of 1.0 with- 
out evaporation. Thus, card 0060 will read: 

0060 -2 0 0.0 150.0 1.0 0 


Four 0061 cards are expected, recording the four y posi- 
tions of the particles and the proper location of the spray boom 
relative to the surface. These cards will read: 


0061 

-3.0 

-0.75 

0061 

0.5 

-0.75 

0061 

2.6 

-0.75 

0061 

4.8 

-0.75 


The AGDISP code can now be run. The terminal output for this 
example is shown in Figure 3 and the first and last pages of the 
printer output are shown in Figure 4. Plotting vertically on a 
Tektronix 4025 with convenient scales to include the particle mean 
trajectories and standard deviations gives the trajectory plot in 
Figure 5 and the deposition plot in Figure 6. Invoking AGLINE for 
the maximum f igure-of-merit value results in the distribution plot 
shown in Figure 7. 
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7. AGDISP FLOW CHART 


The AGDISP, AGPLOT and AGLINE codes consist of mainline 
sections that process user-supplied input data and data file 
information, initializing the program fields and output files. 

From these main routines, subroutines are called to provide 
specific information pertinent to the solution of the particle 
equations. The simplified models for particle time scale, evap- 
oration time scale, particle- turbulence correlation, etc. , are 
placed in separate subroutines for ease of modification. This 
section of the User Manual examines the flow charts to these 
three codes. 

The AGDISP code is configured to run interactively on the 
Langley computer, with a size under 60K octal words. Input comes 
from a card reader file on unit 4, with terminal output written 
on unit 6, and printer output written on unit 9. The code gener- 
ates a trajectory plot file on unit 8, to be used subsequently by 
the AGPLOT code. When an externally- computed fluid- turbulence 
field is supplied by the user, this data file is entered into the 
AGDISP code on unit 10. The unit 8 and unit 10 data files are 
binary. 

The AGPLOT code is also configured to execute interactively, 
with a size under 50K octal words. All input data is entered by 
the user at the terminal on unit 5, and all output is written on 
the terminal in graphical form on unit 6. Only Tektronix and 
Langley computer center graphics are available. The data file 
generated by AGDISP serves as the input data file for AGPLOT; this 
file is entered on unit 8. 

The AGLINE code executes interactively, with a size also 
under 50K octal words. Input data is entered by the user at the 
terminal on unit 5, with trajectory data from the AGDISP plot file 
on unit 8. Terminal output is written on unit 6 and printer out- 
put on unit 9. 

The flow charts for the AGDISP, AGPLOT and AGLINE codes are 
given in Figures 8, 9 and 10, respectively. These charts attempt 
to show the overall details of the respective codes , indicating 
what each subroutine does; and supplementing the discussion of 
Section 8. Procedure files used to execute these programs are 
illustrated in Tables 1, 3 and 4. 
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8. AGDISP SUBROUTINES 


This section of the User Manual examines the code comprising 
each subroutine in the AGDISP, AGPLOT and AGLINE computer programs. 
While Figures 8, 9 and 10 give an overview of the operation of the 
programs, this section is intended to summarize the contents of 
the code. The routines are examined in alphabetical order. 


AGDISP is the mainline program computing particle motion by a 
Lagrangian technique. Most of the details of AGDISP 

are covered earlier in the User Manual, including data 
card entry, verification, initialization and calculation. 
Each data card has its own section in AGDISP wherein 
valid data entry is checked and initialization is per- 
formed. Card order is important, and is maintained by 
a series of error flags L20, L22, L25 , etc. Free-format 
entry demands that all data cards contain the appropriate 
data; otherwise, an error will be trapped. Inadvertently 
reaching the end of the data deck also generates an 
error. All input cards generate output onto the terminal. 
After completing the input, any possible warning condi- 
tions are flagged before selecting the integration step 
size equal to one-half the particle time constant (this 
time step may be smaller if evaporation is included in 
the calculation). Also, the code has been configured 
to adjust the output frequency such that line printer 
output occurs approximately every one-tenth second of 
simulated time. Integration is then performed, after 
which the final particle ground deposition is computed, 
taking evaporation into account, and the input card 
deck is written to the generated plot file for typing 
back to the user when plotting the trajectories and 
deposition with AGPLOT. 

AGBZD evaluates the derivatives for the time-dependent solu- 
tion of the Betz roll up methodology found in Ref. 1. 

The input parameters are the present values of the 
variables, the derivatives (to be computed), the start- 
ing pointer, maximum pointer, end pointer and total entry 
pointer of the rolling up vortex circulation distribution. 
The variables are R , the radius of the vortex core, 
and T , the circulation of the vortex. At t = 0 , 

R is zero and the singularity is treated by taking the 
appropriate limit (see Ref. 1). 

AGBZG initializes the Betz roll up procedure using the user- 

inputted circulation distribution stored in common block 
BETZ. The first section of the subroutine computes the 
spatial derivative dr/dy of the circulation V and 
sets a noise level of 0.5% of peak derivative value. 

The circulation distribution is then reexamined to locate 
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AGBZI 


AGBZQ 


AGBZT 

AGCOR 


where dr/dy is maximum (these become positions where 
vortices start rolling up) and where dr/dy is minimum 
(these become the end positions of the vortex circula- 
tion patterns). The noise level is used here to correct 
for slightly nonuniform data (a frequent occurrence when 
discrete data is entered) . The circulation pointers to 
start, maximum and end for up to four vortices are then 
set (more than four vortices invoke an error exit from 
AGDISP) . The initialization is completed by evaluating 
the derivatives at zero time and initializing all of the 
parameters pointing in the AGDISP code to vortex center 
location, strength, canopy effect and unrolled sheet 
effect. The code then establishes a vortex-dependent 
time step DTV vector based on 1% of the rollup time 
constant x R - 2irR^/r . 


integrates the Betz equations across the time step DELT 
entered from the AGDISP code. The first section of the 
subroutine establishes the time step and number of steps 
to integrate based on the computed values of DTV and 
the vortices not yet rolled up completely. Each vortex 
is treated separately by first predicting its new values 
of R and r , then solving for the derivatives at these 
values and correcting the solution twice. This double 
correction enhances solution accuracy (Ref 11). The 
y centroid of the rolling up vortex is computed based 
upon how much circulation has rolled into the vortex 
by this time. The vertical centroid z is not incre- 
mented. The incremental movement of the vortex for this 
Betz time step is added to the current position of the 
vortex (this position may be influenced by the other 
background features incorporated in AGDISP) , and fully 
rolled-up vortices are flagged with zero sheet strength. 
Lastly, the unrolled sheet lengths are determined. 


is a Gauss -Legendre integration routine used to find the 
area under the inputted discrete function between input- 
ted starting and ending points (Ref. 12). 


is a linear interpolator extracting the value of the in- 
putted discrete function at the inputted position. 

evaluates the particle-turbulence correlations at time 
T , with particle time constant DTAU and turbulence 
time constant WTAU . The equations are a result of the 
assumption of a von Karman spectral distribution as des- 
cribed in Ref. 1. The constants XK1 , XK2 and XK3 
take their limiting values whenever the two time constants 
are within 0.5% of each other. The correlations <ux> 
and <uv> are then evaluated, and these become both 
<uy> , <uz> and <uv> , <uw> respectively. 
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AGDEC 


AGDIF 


AGINT 


AGLQD 


AGLQS 

AGMAT 


computes the decay rate based on the particle time 
constant DTAU , with the effect of evaporation being 
a reduction of the time constant consistent with the 
smaller size particle. When the particle size reaches 
the termination size, the end-of-run flag LENDF is 
set. The decay rate is simply the sum of two terms, be- 
cause of the simple evaporation model built into the 
AGDISP code. The drag law of Langmuir and Blodgett 
(Ref. 4) is used. 


evaluates the derivatives of the Lagrangian particle 
equations discussed in detail in Ref. 1. The time con- 
stant is evaluated, and the equations for each particle 
in the simulation are evaluated. The first and sixth 
elements of the vector XV are the particle position 
(y,z) . At this point, the mean velocity components V 
and W are determined, the background turbulent time 
constant is determined and the analytic particle turbu- 
lence correlations are computed. The equation deriva- 
tives are formed, with gravity equal to 9.8 m/sec . 

monitors the integration of the particle equations. 
Initially, all particle positions are stored and initial 
derivatives determined. Integration then proceeds step 
by step to the maximum inputted time. The WAKE plot file, 
Betz, canopy, propeller and helicopter are updated where 
applicable, particle positions are predicted, derivatives 
evaluated, and the positions are twice corrected. Part- 
icle solutions terminate at the ground. Since the vort- 
ices move under the influence of one another and their 
ground images, their positions are also adjusted. Any 
particles that have impacted the surface are flagged and 
their final sizes are computed. Termination is checked 
for evaporation, maximum time, and ground impact; and 
plot save is invoked. If termination exists, transfer is 
returned to AGDISP; otherwise another time step is taken. 


is a linear decomposition routine transposed from the 
IMSL scientific subroutine package (Ref. 13) , used to 
decompose the six-by-six matrix in the superequilibrium 
methodology. 

is the linear substitution routine from Ref. 13, used to 
evaluate the six super equilibrium unknowns. 


fills the six-by-six matrix array for the unknowns 
<uu> , <w> , <ww> , <uv> , <uw> and <vw> by super- 
equilibrium (Ref. 9 ) for a given value of q^ . The 
linear equations are solved by AGLQD and AGLQS to 


determine the difference between 
<uu> 4- <w> + <ww> . 


,2 


and 
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AGPAC 


AGP AD 


AGPRP 


AGRTF 


AGSAV 


AGS UP 


computes the vortex circulation reduction resulting from 
interactions with a canopy. Essentially, the scrubbing 
of the vortex acts as a drag on the wake flow field. 

This drag translates into an effective vortex strength 
smaller than the noncanopy value. Reference 1 presents 
the details. In this subroutine, vortices are checked 
for whether they have entered the canopy; if so, an 
integration across the portion of the vortex inter- 
acting with the plant area density is taken, a simple 
time integration is performed, and the vortex strength 
reduction factor is computed. The reduction factors 
FACR and FACL modify the vortex strengths, unrolled 
sheet strengths, propeller swirl and helicopter effects. 

computes the displacement thickness of the plant canopy 
and initializes the effective circulation integrals. 

The displacement thickness is assumed to represent the 
effective surface roughness the canopy gives the atmos- 
pheric flow above it. Within the canopy, the crosswind 
velocity and turbulence are assumed to behave linearly 
with height, an approximation consistent with the nu- 
merical prediction of turbulent behavior within plant 
canopies (Ref. 8). 

updates the turbulence, swirl velocity and wake radius 
behind the propeller, using the approximations discus - 
s ed in Ref . 1 . 

is a modified root- finder subroutine taken from the 
IMSL scientific subroutine library (Ref. 13). Essen- 
tially, the routine performs higher-order interpolation 
to yield a solution for q 2 which agrees accurately 
with <uu> + <w> + <ww> obtained by solving the six- 
by-six matrix set. The solution speed is enhanced by 
starting the iteration at a value of q2 close to the 
answer . 

writes the step integration results from AGDISP to the 
plot file (the variables y , z , <yy>^ , <zz>^ ) and 
the line printer (all variables for all nonimpacted 
particles) . 

is the driving routine for superequilibrium. Given the 
four spatial derivatives in the AGDISP code (dU/dy 
and dU/dz are zero) and the local scale length A , 
the superequilibrium equations are iterated for the 
values of <uu> , <w> and <ww> . The maximum grad- 
ient is used to normalize the solver, an approximate 
result is encountered, and then an accurate solution 
is found by stepping ih q until a zero crossing is 
found by invoking the root finder AGRTF to obtain 
a value of q accurate to four digits or a normalized 
value of q2 - <uu> - <w> - <ww> within ± 0.001 m^/sec^ 
of zero. This result produces the values of turbulent 
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AGSVE 


energy unnormalized by A and the normalizing velocity 
gradient. Any (rare) error trap returns laminar results 
in an effort to reduce the impact of the background 
solution on the particle equations. 

computes the incremental background velocity arising from 
the unrolled Betz sheets. In this case, the sheets are 
represented by a constant circulation sheet of vorticity, 
and the resultant flow is analytically determined. The 
singularity in the vicinity of the sheet is controlled 
by keeping all data evaluation at a distance of 1% of 
the sheet length with linear falloff within 10%. As the 
Betz procedure continues , the sheet becomes shorter un- 
til it is finally rolled up. 

AGTUR determines the turbulent scale length A and turbulent 

components q z , <w> and <ww> at the position (y,z). 
The scale is essentially the smaller of 0.65z, the maxi- 
mum inputted scale A max , or six-tenths the distance to 

the nearest vortex. The turbulence comes from super- 
equilibrium or the fixed inputted value of q^ ax (modi- 
fied by the propeller wake) . The values will be modified 
if the particle position is within the canopy (a linear 
correction to the height of the canopy). Alternately, 
the scale and correlations could in all or part come 
from a WAKE plot file. 

AGVCH examines each vortex location and propeller and helicopter 
centerlines and corrects for their movement in the given 
time step DELT 

AGVEL computes the mean velocity components V and W at the 
position (y , z) . Each vortex, its reflecting image 
across y = 0 and their images across z = 0 are used 
to compute the overall velocity increment. The standard 
potential vortex velocity field r / 2 ttt , producing a 
velocity normal to the radius vector, is broken into its 
(y,z) components and modified by FACR and FACL for the 
presence of a canopy. The unrolled sheet effect (with 
image sheets) is added, as are the helicopter, propeller 
and mean crosswind modified by the canopy. Alternately, 
the WAKE plot file is quizzed to return the appropriate 
velocity components. 

AGWKI linearly interpolates the data array read from the WAKE 
plot file for the variable desired (V , W , q z , A , 

<vv> or <ww>) at the position (y,z) . A warning mes- 
sage is written the first time spatial extrapolation 
beyond the WAKE plot file grid coordinate is invoked, 
even though extrapolation may continue indefinitely 
thereafter . 
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AGWKR 


AGWKS 


AGPLOT 


AGCHK 


AGDRW 


AGGRD 


AGMDL 


recovers the WAKE plot file data for the desired time T . 
Interpolation between time steps is performed every one- 
tenth of the interval; otherwise the data arrays in 
common block WAKE are not updated. The file is rewound, 
the initial data skipped, and the pertinent time step 
bracketed. Linear interpolation is then performed, and 
the next time check is computed. When the end of the 
WAKE plot file is reached, the last time entry data is 
used for the duration of the AGDISP run. 

processes the WAKE plot file header information. The 
file is binary and is read using sequential reads. The 
plot file comment of nineteen words is followed by the 
number of variables, a list of the variable names, the 
number of y data points, the values of the y mesh, 
the number of z data points and the values of the 
z mesh. Thereafter the data consists of a time entry 
followed by all y data values for each z mesh row 
for each variable. The structure of the WAKE plot file 
used in test case #5 is shown in Figure 11, where the 
binary information has been printed. Details are given 
in Section 3. 

is the mainline program for plotting the results from 
AGDISP , quizzing the interactive user about the details 
of the desired plots and plotting the trajectory pattern 
and ground deposition. 

verifies the appropriateness of the scaling information 
requested by AGPLOT of the interactive user. Thus, 
the scale must be monotonically increasing with a posi- 
tive spacing between scale marks and an integer number 
of scale marks, not to exceed ten. 

translates the data to be plotted into the graphical 
grid point coordinates representing the data, given the 
mesh sizes to be plotted on the screen. This transformed 
data is plotted by a call to AGMDL . 

computes the ground deposition pattern, quizzes the user 
for plot scale information, and plots the resulting dis- 
tribution. Each impacting particle in the AGDISP simu- 
lation hits the surface at a known y position with a 
known variance <yy> . Adding all particle locations 
with their known impact sizes produces the deposition 
pattern along the surface. The scales are established 
and the curve is plotted. 

plots the graphical screen data points, either directly 
for Tektronix 4025 or with CDC plotting routines for 
Tektronix 4010 and 4014. 
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AGPLT is the plotting subroutine that uses either Tektronix 
4025 or CDC graphics calls to perform the actual plot- 
ting of trajectory data. The screen is first cleared; 
then axes are established and plotted and a wing sche- 
matic is drawn if appropriate. The plot file is read 
so that nine points are plotted at a time, drawing 
the mean particle paths (including the image if desired) 
the standard deviation paths and tag markers at 
the desired intervals. At the conclusion of the plot, 
the screen is copied. 

AGSPD uses the time dependent values of y , z , <yy>% and 
<zz>^ to determine the plotting coordinates of the 
particle standard deviation. The deviation is found 
normal to the instantaneous particle direction, and the 
coordinates are adjusted accordingly. 

Simple FORTRAN WRITES are needed to cue the Tektronix 4025 
to construct a graphics area, write into it and produce a printed 
copy of the screen. The default character is the exclamation 
mark (!) and the useful graphics commands are the following: 

BEL - rings the terminal bell (signals end-of-plot) . 

GRA - establishes a graphics area. 

HCO - invokes a copy of the screen. 

LIN - sets the line segment character (1 is solid, 3 is dashed). 

MON - transfers command response to the terminal monitor space. 
RDO - positions the workspace. 

VEC - draws line segments in the workspace. 

WOR - establishes a terminal workspace and transfers command 
response to it. 

For the Tektronix 4010 and 4014 subroutine calls must be made 
using the CDC-Langley software to accomplish the graphics. The 
subroutines invoked by AGPLOT are the following: 

BELL - rings the terminal bell. 

DRASE - erases the screen, returns to character mode and sets 

the cursor to the home position (upper left hand corner). 

DRWABS - draws a solid line from the current cursor position to 
the screen coordinate supplied. 

DSHABS - draws a dashed line (with defined characteristics) 

from the current cursor position to the screen coordin- 
ate supplied. 
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HDCOPY 


MOVABS 

PSEUDO 

TSEND 

AGLINE 


AGCHK 

AGDRW 

AGELP 

AGEQD 


AGMDL 

AGPLT 


invokes a copy of the screen. 

moves the cursor from wherever it is to the screen 
coordinate supplied (no line is drawn) . 

erases the screen, establishes a plotting buffer, 
enters graph mode and positions the cursor at (0,0). 

flushes the plotting buffer, thereby finishing the 
current plot. 

is the mainline program for interpreting the AGD1SP 
plot file to construct an equivalent Gaussian profile 
distribution from the multiple-particle solution, 
compatible with the FSCBG code (Ref. 10). Three 
passes are made through the plot file. The first 
locates the AGDISP input data cards (at the end of 
the file) and displays them on the screen. The 
second pass computes the equivalent Gaussian profile 
distribution and figure-of -merit given the multiple- 
particle solution at every saved time in the plot 
file. The third pass locates the desired sequence 
number or figure-of -merit and plots the final 
distributions . 

verifies the appropriateness of the scaling informa- 
tion for AGLINE as it does for AGPLOT. 

translates the data to be plotted into grid coordin- 
ates as for AGPLOT. 

sets up the plots of particle isopleths at one-eighth, 
one- fourth, one-half, one and twice the standard 
deviation. 

computes the equivalent Gaussian distribution of the 
multiple-particle solution (including the particle 
solution on the other half of a half-plane solution) . 
The mean horizontal and vertical positions and stan- 
dard deviations are summed over all particles in the 
solution plane. 

plots the graphical screen data as for AGPLOT. 

is the plotting subroutine, working similar to the 
AGPLT subroutine in AGPLOT, plotting the particle and 
equivalent distributions and copying the screen. 
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9. AGDISP TEST CASES 


Six test cases are included in the User Manual to illu- 
strate the variety of data entry available with AGDISP. These 
examples are included to demonstrate the structure of the input 
card deck and the data entry requirements of each type of data 
card and option. They are not meant to suggest that certain op- 
tions go with certain other options, nor are all options covered. 
For example, only one WAKE plot file run is shown, in this case 
for a plot file containing V , W and with A and <w> 

and <ww> determined algebraically. Extension to more complete 
files has been demonstrated but is not included. 

The six examples are meant to illustrate the following fea- 
tures : 

1. A rectangularly loaded wing, with fixed background 
turbulence, full plane results with code-determined 
particle positions, including a centerline data point 
and a propeller. 

2. A rectangularly loaded biplane wing, with superequili- 
brium invoked, half -plane results with code-determined 
particle positions, no centerline data point, with 
evaporation. 

3. Circulation loading with Betz, a canopy, half-plane 
with fixed turbulence and the particles user- inputted . 

4. A single particle release into a crosswind, fixed 
turbulence in a full plane. 

5. A WAKE plot file background condition for particles 
user-inputted in a full-plane configuration. 

6. A helicopter input for a fixed turbulence background 
with code-determined particle positions in a half-plane 
simulation. 

The simulation time is kept short since the essential fea- 
tures of the input data decks and the general character of the 
evolving solution are what are sought here. Each test case makes 
use of the documentation in Section 3 to construct the 
input decks. Figures 12-17 illustrate the six test cases. Each 

figure is a reproduction of the terminal output (showing the 
reproduction of the input deck and the AGDISP echo responses) 
and the first and last output pages (showing the final time step 
for the nonimpacted particles and giving the deposition informa- 
tion) . An example of plotted particle and standard deviation 
results and the corresponding ground deposition plot may be found 
for the example test case in Section 6. 
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TABLE 1 


Typical AGDISP Interactive Procedure File 


Command Line 


Description 


MAP, OFF. 

AS SIGN, MS, OUTPUT. 
ASSIGN,TT,TAPE6. 

GET,TAPE4=INPUT0 . 


GET, TAPE 10=WAKE0. 


GET, OBDISP . 

OBDISP. 

REWIND, TAPE8,TAPE9. 
REPLACE .TAPE 8=PL0T0 . 


REPLACE , TAPE 9=0UTPUT0 . 


Linkage editor load map output is disabled. 

Program underflows are written to the file 
OUTPUT. By assigning OUTPUT to mass storage, 
these warning messages are suppressed. 

AGDISP writes all interactive output to unit 
number 6. This unit is assigned to the time- 
sharing terminal and the output appears on the 
screen. 

AGDISP reads all input data cards from unit 
number 4. This unit is equated to the default 
input file name INPUTO and is retrieved from the 
permanent directory as the local file TAPE4. 

AGDISP reads wake plot file input data from unit 
number 10. This unit is equated to the default 
wake file name WAKEO and is retrieved from the 
permanent directory as the local file TAPE10. 
NOTE: This GET is issued only when the wake 
plot file code option is invoked. 

The AGDISP object module OBDISP is retrieved 
from the permanent directory. 

Load and execute the AGDISP object module. 

Output written by AGDISP to unit numbers 8 
and 9 are positioned at their initial records. 

AGDISP writes trajectory plot file data results 
to unit number 8 . This unit is equated to the 
default plot file name PL0T0 and transferred 
from the local directory to the permanent 
directory, overwriting existing data. 

AGDISP writes printer output data results to 
unit number 9 . This unit is equated to the 
default output file name 0UTPUT0 and 
transferred from the local directory to the 
permanent directory, overwriting existing data. 
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TABLE 2 


0000 

0010 

0000 

0015 

0000 

0020 

0000 

0021 

0000 

0022 

0000 

0025 

0025 

0025 

0025 

0025 

0000 

0028 

0000 

0030 

0000 

0040 

0000 

0050 

0000 

0055 

0055 

0055 

0055 

0000 

0060 

0000 

0061 

0061 

0061 

0061 

0000 

0062 

0000 

0065 


Example Card Input to AGDISP 


card 

0010 

is 

the program data card 


10.0 

card 

1 

0015 

is 

the terrain slope card 


0.0 

card 

0020 

is 

the aircraft characteristics 

data card 

2 0 

6.5 

4.0 

35.0 1 


card 

0021 

is 

the biplane data card 


2.0 

card 

0.8 0 
0022 

.6 

specifies vortex circulation 


25.0 

card 

0025 

specifies Betz wing circulation 

data 


0.0 10.0 
2.0 10.0 

3.0 20.0 

4.0 15.0 

-5.0 0.0 

card 0028 is the cross wind velocity data card 

3.0 3.0 0.01 

card 0030 is the helicopter data card 

15000.0 0.05 

card 0040 is the propeller data card 
0.1 16.0 0.8 500.0 2.0 0.75 

card 0050 is the turbulence data card 
0 1.0 30.0 


card 

0055 

specifies the canopy 

profile data 

0.0 

0.05 



0.3 

0.1 



0.7 

0.3 



-1.5 

0.0 



card 

0060 

is the particle data 

card 

-3 -1 

0.0 

100.0 1.0 1 



card 0061 enters discrete particle locations 
1.0 0.0 
2.5 0.0 

4.0 0.0 

0.0 0.0 

card 0062 provides particle initial conditions 
0.0 -0.5 0.5 0.0 

card 0065 is the evaporation data card 

5.0 50.0 
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TABLE 3 


Typical AGPLOT Interactive Procedure File 


Command Line 


Description 


MAP, OFF. 

ASSIGN, MS, OUTPUT. 
ASSIGN,TT,TAPE6. 

GET , TAPE8=PLOTO . 


Linkage editor load map output is disabled. 

Program underflows are written to the file 
OUTPUT. By assigning OUTPUT to mass storage, 
these warning messages are suppressed. 

AGPLOT writes all interactive output to unit 
number 6. This unit is assigned to the time- 
sharing terminal and the output appears on the 
screen. 

AGPLOT reads all trajectory plot data from 
unit number 8 . This unit is equated to the 
default trajectory plot file name PLOTO and 
is retrieved from the permanent directory as 
the local file TAPE8. 


ATTACH ( LI BFTEK/UN=LIBRARY) To write graphics data to the Tektronix 401x, 

basic system plotting subroutines called by 
AGPLOT must be made available to the linkage 
editor. These routines are located within the 
system library LIBFTEK, and this library is 
made a local library. 

GET,OBPLOT. The AGPLOT object module OBPLOT is retrieved 

from the permanent directory. 

LDSET(LIB=LIBFTEK) The linkage editor is instructed to search 

LIBFTEK for the graphics subroutines (BELL, 
DRASE, DRWABS, DSHABS , HDCOPY, MOVABS, 

PSEUDO, and TSEND) . 

OBPLOT. Load and execute the AGPLOT object module. 
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TABLE 4 


Typical AGLINE Interactive Procedure File 


Command Line 


Description 


MAP, OFF. 

AS SIGN, MS, OUTPUT. 


ASSIGN, TT,TAPE6. 


Linkage editor load map output is disabled. 

Program underflows are written to the file 
OUTPUT. By assigning OUTPUT to mass storage, 
these warning messages are suppressed. 

AGLINE writes all interactive output to unit 
number 6. This unit is assigned to the time- 
sharing terminal and the output appears on the 
screen. 


GET, TAPE 8=PL0T0 . AGLINE reads all trajectory plot data from 

unit number 8 . This unit is equated to the 
default trajectory plot file name PLOTO and 
is retrieved from the permanent directory as 
the local file TAPE8. 


ATTACH(LIBFTEK/UN=LIBRARY) To write graphics data to the Tektronix 401x, 

basic system plotting subroutines called by 
AGLINE must be made available to the linkage 
editor. These routines are located within the 
system library LIBFTEK, and this library is 
made a local library. 

GET,OBLINE. The AGLINE object module OBLINE is retrieved 

from the permanent directory. 

LDSET(LIB=LIBFTEK) The linkage editor is instructed to search 

LIBFTEK for the graphics subroutines (BELL, 
DRASE, DRWABS , DSHABS , HDCOPY, MOVABS, 

PSEUDO, and TSEND) . 


OBLINE. Load and execute the AGLINE object module. 

REWIND, TAPE9 . Output written by AGLINE to unit number 9 

is positioned at its initial record. 

REPLACE , TAPE 9=0UTPUT0 . AGLINE writes printer output data results to 

unit number 9. This unit is equated to the 
default output file name OUTPUTO and 
transferred from the local directory to the 
permanent directory, overwriting existing data. 
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Figure 2. Airplane configuration for the example test case of Section 


NASA AGDISP (MOD 2.0) PROGRAM RESULTS 


1: 

0000 

EXAMPLE TEST CASE 

2: 

0010 

10.0 2 

3: 

0020 

2 1 6.0 3.75 35.0 

4: 

0022 

33.6 

5: 

0028 

2.0 5.0 0.02 

6: 

0040 

0.1 18.0 0.8 500.0 

7: 

0050 

0 0.0 30.0 

8: 

0060 

-2 0 0.0 150.0 l.i 

9: 

0061 

-3.0 -0.75 

10: 

0061 

0.5 -0.75 

11: 

0061 

2.6 -0.75 

12: 

0061 

4.8 -0.75 


FULL-PLANE CALCULATION 

AIRCRAFT SEMI -SPAN/DISK RADIUS: 0.60000E+01 M 

FLIGHT SPEED: 0.35000E+02 M/SEC 

NOMINAL RELEASE HEIGHT: 0.37500E+01 M 

RECTANGULARLY LOADED WING WITH GAMMA: 0.33600E+02 M' 


: 2/SEC 


CROSS-WIND VELOCITY 
Z 
ZO 

PROPELLER HUB HEIGHT 
RADIUS 
SWIRL VELOCITY 
TURBULENCE 


0 . 20000E+01 M/SEC 
0 . 50000E+01 M 
0 . 20000E-01 M 
0 . 25000E+00 M 
0 . 15000E+01 M 
0 . 46835E+01 M/SEC 
0.31810E+01 (M/SEC)**2 
TURBULENCE FIXED VALUE: 0.11071E+00 (M/SEC)**2 

SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

TOTAL NUMBER OF PARTICLES: 4 

DIAMETER: 0.15000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

INITIAL TIME STEP: 0.35100E-01 SEC 

MAXIMUM TIME: 0.10000E+02 SEC 


TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 

TIME 


Figure 3. 


0.3866E+00 

SEC 




0 . 8143E+00 

SEC 




0 . 1261E+01 

SEC 




0. 1713E+01 

SEC 




0 . 2168E+01 

SEC 

TIME 

0 . 6983E+01 

SEC 

0 . 2627E+01 

SEC 

TIME 

0 . 7312E+01 

SEC 

0 . 3032E+01 

SEC 

TIME 

0 . 7641E+01 

SEC 

0 . 3362E+01 

SEC 

TIME 

0. 7971E+01 

SEC 

0 . 3691E+01 

SEC 

TIME 

0.8300E+01 

SEC 

0 . 4020E+01 

SEC 

TIME 

0.8629E+01 

SEC 

0.4349E+01 

SEC 

TIME 

0 . 8958E+01 

SEC 

0 . 4679E+01 

SEC 

TIME 

0 . 9287E+01 

SEC 

0 . 5008E+01 

SEC 

TIME 

0 . 9617E+01 

SEC 

0.5337E+01 

SEC 

TIME 

0 . 9946E+01 

SEC 

0 . 5666E+01 

SEC 




0 . 5995E+01 

SEC 

INTEGRATION COMPLETE 

0 . 6325E+01 

SEC 




0 . 6654E+01 

SEC 

NASA 

AGDISP (MOD 2 

.0) 


Terminal output from the AGDISP run of the Example 
Test Case constructed in Section 6. 
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Figure 5. Trajectory plot of the particle means and standard deviations 
for the Example Test Case of Section 6. 
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Figure 6. Ground deposition pattern for the Example Test Case of Secti 
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Figure 7. Equivalent Gaussian distribution at maximum f igure-of -merit 
for the Example Test Case of Section 6. 
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Figure 8. AGDISP program flow chart 
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Figure 8. (Cont’d) 
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Figure 8. (Cont’d) 
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Figure 9. AGPLOT program flow chart 


47 


II 



AGLINE 


! START | 

L ! 


READ ADGISP DISK FILE TO LOCATE CASE CARDS 


I 

j 


READ AGDISP DISK FILE TO COMPUTE EQUIVALENT 
GAUSSIAN DISTRIBUTION AND FIGURE-OF-MERIT 
FOR ALL ENTRIES IN DISK FILE 


INPUT: SEQUENCE NUMBER OR FIGURE-OF-MERIT? | 

... i 


| READ AGDISP DISK FILE TO FIND DESIRED POSITION | 

i ... ...... j 


| INPUT: y SCALE DATA? ] AGCHK 

t INPUT: z SCALE DATA? | AGCHK 

i i 


PLOT DISTRIBUTIONS: j AGPLT 

INITIALIZE SCREEN S SCALES | 

ENTER AIRCRAFT SCHEMATIC | AGELP 

PLOT PARTICLE PROFILES (DASHED) | AGDRW 

PLOT EQUIVALENCE PROFILES (SOLID) j AGMDL 

— — i 


| HARD COPY j 

i i 


i 1 

j END j 

I i 


Figure 10. AGLINE program flow chart 
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Record Contents 


Description 


TEST WITH VELOCITIES AND QQ WAKE plot file header record - 76 characters long. 

3 Number of variables in plot file - the maximum 

number is six: horizontal and vertical 
velocities V and W in m/sec; turbulence QQ in 
(m/sec)**2; turbulent macroscale SL in m; and 
horizontal and vertical turbulence components 
W and WW in (m/sec)**2. 

V W QQ Descriptive titles of the variables - three in 

this case - 4 characters per title. Verification 
of number of variables in plot file with flags 
in the AGDISP input data deck is made. 

3 Number of Y (horizontal) data points in plot 

file - 16 is the maximum permitted. 

0.0 20.0 40.0 The three Y position values in meters. If there 

are more Y values, they go here also. 

3 Number of Z (vertical) data points in plot 

file - 16 is the maximum permitted. 

0.0 20.0 40.0 The three Z position values in meters. If there 

are more Z values, they go here also. 

0.0 The first time in the plot file: 0.0 seconds. 

0.0 20.0 40.0 The values of the first variable (V) for the 

first Z position (0.0) for all Y positions 
in the file. Additional values of V go here 
if there are more Y positions. 


0.0 

20.0 

40.0 

The values of V for the second Z position (20.0) 
for all Y positions in the file. 

0.0 

20.0 

40.0 

The values of V for the third Z position (40.0) 
for all Y positions in the file. Additional 
values of V for additional Z positions go here 
if there are more Z positions. 

0.0 

0.0 

0.0 

The values of the second variable (W) for the 
first Z position for all Y positions. 

-20.0 

i 

to 

o 

o 

i 

to 

o 

o 

The values of W for the second Z position for 
all Y positions. 

o 

0 

1 

o 

0 
vf 

1 

o 

0 

1 

The values of W for the third Z position for 


all Y positions. Additional values of W go 
here for additional Z positions. 


Figure 11. Structure of the WAKE plot file used in test case #5 
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1.0 

1.0 

1.0 

The 

the 

values of the third variable (QQ) for 
first Z position for all Y positions. 

1.0 

1.0 

1.0 

The 

all 

values of QQ for the second Z position for 
Y positions. 

1.0 

1.0 

1.0 

The 

values of QQ for the third Z position for 


all Y positions. Additional values of QQ go 
here for additional Z positions. 


Additional variable data (for SL, W and WW) 
follow the same pattern and go here if 
applicable . 


2.0 



The second time in the plot 

file: 2.0 seconds. 

0.0 

10.0 

20.0 

The V data corresponding to 

the second time. 

0.0 

10.0 

20.0 



0.0 

10.0 

20.0 



0.0 

0.0 

0.0 

The W data corresponding to 

the second time. 

-10.0 

-10.0 

-10.0 



-20.0 

-20.0 

-20.0 



0.05 

0.05 

0.05 

The QQ data corresponding to the second time. 

0.05 

0.05 

0.05 



0.05 

0.05 

0.05 






Additional variable data go 

here if applicable 

5.0 



The third time in the plot : 

file: 5.0 seconds. 

0.0 

4.0 

8.0 

The V data corresponding to 

the third time. 

0.0 

4.0 

8.0 



0.0 

4.0 

8.0 



0.0 

0.0 

0.0 

The W data corresponding to 

the third time. 

-4.0 

-4.0 

-4.0 



-8.0 

-8.0 

-8.0 



0.02 

0.02 

0.02 

The QQ data corresponding to the third time. 

0.02 

0.02 

0.02 



0.02 

0.02 

0.02 






Additional variable data go 

here if applicable 


Additional times and variable data go here 
if applicable, repeating the pattern 
established above. 


-1.0 The WAKE plot file termination signal, a 

negative time. 

Figure 11. (Cont'd) 
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NASA AGDISP (MOD 2.0) PROGRAM RESULTS 


1: 

0000 

AGDISP TEST CASE # 1 


2: 

0010 

10.0 2 




3: 

0020 

2 0 6.5 

4.0 

35.0 

0 

4: 

0022 

25.0 




5: 

0040 

0.1 16.0 

0.8 

500.0 

2 

6: 

0050 

0 1.0 30. 

0 



7: 

0060 

3 1 0.0 

100. 

o 

rH 

O 

0 


FULL-PLANE CALCULATION 

AIRCRAFT SEMI -SPAN/DISK RADIUS: 0.65000E+01 M 

FLIGHT SPEED: 0.35000E+02 M/SEC 

NOMINAL RELEASE HEIGHT: 0.40000E+01 M 

RECTANGULARLY LOADED WING WITH GAMMA: 0.25000E+02 M**2/SEC 

PROPELLER HUB HEIGHT: 0.75000E+00 M 

RADIUS: 0 . 20000E+01 M 

SWIRL VELOCITY: 0.18060E+01 M/SEC 

TURBULENCE: 0.84092E+00 (M/SEC )**2 

TURBULENCE FIXED VALUE: 0.10000E+01 (M/SEC)**2 

SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

TOTAL NUMBER OF PARTICLES: 7 

DIAMETER: 0.10000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

INITIAL TIME STEP: 0.15600E-01 SEC 

MAXIMUM TIME: 0.10000E+02 SEC 

TIME: 0 . 7066E+00 SEC 

TIME: 0.1422E+01 SEC 

TIME: 0 . 2148E+01 SEC 

TIME: 0 . 2881E+01 SEC 

TIME: 0 . 3615E+01 SEC 

TIME: 0 . 4349E+01 SEC 

TIME: 0 . 5084E+01 SEC 

TIME: 0 . 5820E+01 SEC 

TIME: 0 . 6491E+01 SEC 

TIME: 0 . 7024E+01 SEC 

TIME: 0 . 7557E+01 SEC 

TIME: 0 . 8089E+01 SEC 

TIME: 0 . 8622E+01 SEC 

TIME: 0 . 9154E+01 SEC 

TIME: 0 . 9687E+01 SEC 

INTEGRATION COMPLETE 

NASA AGDISP (MOD 2.0) PROGRAM END 


Figure 12a. AGDISP test case #1 terminal output 
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NASA AGDISP (MOD 2.0) PROGRAM RESULTS 
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WW 
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3219E+00 
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3250E+00 
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WW 
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NASA AGDISP (MOD 2.0) PROGRAM RESULTS 


1: 

0000 

AGDISP TEST CASE # 2 

2: 

0010 

2.0 1 

3: 

0020 

2 05.5 4.0 25.0 

4: 

0021 

1.5 1.0 1.0 

5: 

0022 

20.0 

6: 

0050 

-1 0.0 30.0 

7: 

0060 

4 0 0.0 100.0 1. 

8: 

0065 

5.0 50.0 

HALF- 

-PLANE 

CALCULATION 


AIRCRAFT SEMI -SPAN/DISK RADIUS: 0.55000E+01 M 

FLIGHT SPEED: 0.25000E+02 M/SEC 

NOMINAL RELEASE HEIGHT: 0.40000E+01 M 

BIPLANE INCREMENTAL HEIGHT: 0.15000E+01 M 

SPAN FRACTION: 0.10000E+01 

GAMMA FRACTION: 0.10000E+01 

RECTANGULARLY LOADED WING WITH GAMMA: 0.20000E+02 M**2/SEC 

TURBULENCE FROM SUPEREQUILIBRIUM 

SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

TOTAL NUMBER OF PARTICLES: 4 

DIAMETER: 0.10000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

EVAPORATION TEMPERATURE: 0.50000E+01 DEG C 

CUTOFF DIAMETER: 0.50000E+02 MICRONS 

$$$ WARNING: SUPEREQUILIBRIUM TURBULENCE INVOKED 
$$$ WARNING: EVAPORATION INVOKED 
INITIAL TIME STEP: 0.15600E-01 SEC 

MAXIMUM TIME: 0.20000E+01 SEC 

TIME: 0 . 6592E+00 SEC 

TIME: 0 . 1344E+01 SEC 

INTEGRATION COMPLETE 

NASA AGDISP (MOD 2.0) PROGRAM END 


Figure 13a. AGDISP test case #2 terminal output 
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NASA AGDISP (MOD 2.0) PROGRAM RESULTS 

1: 0000 AGDISP TEST CASE # 3 
2: 0010 10.0 1 

3: 0020 0 0 6.16 3.0 30.0 0 

4: 0025 0.0 8.764 

5: 0025 1.733 8.764 

6: 0025 1.959 10.149 

7: 0025 2.185 11.487 

8: 0025 2.412 12.723 

9: 0025 2.638 13.848 

10: 0025 2.865 14.870 

11: 0025 3.091 15.827 

12: 0025 3.318 16.719 

13: 0025 3.544 17.519 

14: 0025 3.770 18.216 

15: 0025 3.997 18.968 

16: 0025 4.223 19.740 

17: 0025 4.450 20.511 

18: 0025 4.676 21.180 

19: 0025 4.903 21.664 

20: 0025 5.129 21.859 

21: 0025 5.356 21.357 

22: 0025 5.582 19.591 

23: 0025 5.809 16.171 

24: 0025 6.035 10.827 

25: 0025 -6.160 0.0 

26: 0050 0 1.0 30.0 

27: 0055 0.0 0.05 

28: 0055 0.3 0.1 

29: 0055 0.7 0.3 

30: 0055 -1.5 0.0 

31: 0060 -4 0 0.0 100.0 1.0 0 

32: 0061 1.0 0.0 

33: 0061 2.5 0.0 

34: 0061 4.0 0.0 

35: 0061 5.0 0.0 

HALF -PLANE CALCULATION 
AIRCRAFT SEMI -SPAN/DISK RADIUS: 0.61600E+01 M 

FLIGHT SPEED: 0 . 30000E+02 M/SEC 

NOMINAL RELEASE HEIGHT: 0.30000E+01 M 

BETZ ROLL UP INVOKED, MAXIMUM ENTRIES: 22 

TURBULENCE FIXED VALUE: 0.10000E+01 (M/SEC)**2 

SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

CANOPY INVOKED, MAXIMUM ENTRIES: 4 

PLANT AREA DENSITY HEIGHT: 0.15000E+01 M 

DISPLACEMENT THICKNESS: 0.73031E+00 M 

TOTAL NUMBER OF PARTICLES: 4 

DIAMETER: 0.10000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

INITIAL TIME STEP: 0.15600E-01 SEC 

MAXIMUM TIME: 0.10000E+02 SEC 

Figure 14a. AGDISP test case #3 terminal output 


TIME: 0 . 6315E+00 SEC 

TIME: 0 . 1324E+01 SEC 

TIME: 0 . 2055E+01 SEC 

TIME: 0 . 2786E+01 SEC 

TIME: 0 . 3487E+01 SEC 

TIME: 0 . 4209E+01 SEC 

TIME: 0 . 4940E+01 SEC 

TIME: 0 . 5666E+01 SEC 

TIME: 0 . 6398E+01 SEC 

TIME: 0 . 7 132E+01 SEC 

TIME: 0.7856E+01 SEC 

TIME: 0.8574E+01 SEC 

TIME: 0 . 9294E+01 SEC 

INTEGRATION COMPLETE 

NASA AGDISP (MOD 2.0) PROGRAM END 
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NASA AGD1SP (MOD 2.0) PROGRAM RESULTS 
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Figure 14b. AGDISP test case #3, first full page of printer output 
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NASA AGDISP (MOD 2.0) PROGRAM RESULTS 


1: 

0000 

AGDISP TEST CASE 

# 4 

2: 

0010 

10.0 2 


3: 

0020 

-2 1 0.0 2.0 

0.0 

4: 

0028 

3.0 3.0 0.01 


5: 

0050 

0 0.0 30.0 


6: 

0060 

0 0 0.0 100.0 

l.l 


FULL -PLANE CALCULATION 
NOMINAL RELEASE HEIGHT: 0.20000E+01 M 

CROSS-WIND VELOCITY: 0.30000E+01 M/SEC 

Z: 0 . 30000E+01 M 

ZO: 0 . 10000E-01 M 

TURBULENCE FIXED VALUE: 0.23349E+00 (M/SEC)**2 

SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

TOTAL NUMBER OF PARTICLES: 1 

DIAMETER: 0.10000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

INITIAL TIME STEP: 0.15600E-01 SEC 

MAXIMUM TIME: 0.10000E+02 SEC 

TIME: 0 . 7204E+00 SEC 

TIME: 0 . 1457E+01 SEC 

TIME: 0 . 2193E+01 SEC 

TIME: 0 . 2929E+01 SEC 

TIME: 0 . 3666E+01 SEC 

TIME: 0 . 4402E+01 SEC 

TIME: 0 . 5 138E+01 SEC 

TIME: 0 . 5875E+01 SEC 

TIME: 0 . 6611E+01 SEC 

TIME: 0 . 7347E+01 SEC 

TIME: 0 . 8083E+01 SEC 

INTEGRATION COMPLETE 

NASA AGDISP (MOD 2.0) PROGRAM END 


Figure 15a. AGDISP test case #4 terminal output 
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0000 AGDISP TEST CASE ff 4 
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NASA 

AGDISP (MOD 

2.0) PROGRAM RESULTS 

1: 

0000 

AGDISP 

TEST CASE # 5 

2: 

0010 

10.0 

2 

3: 

0020 

-1 0 

6.0 3.0 30.0 0 

4: 

0050 

1 0.0 

30.0 

5: 

0060 

-3 -1 

0.0 100.0 1.0 0 

6: 

0061 

-5.0 

0.0 

7: 

0061 

-2.0 

0.0 

8: 

0061 

-1.0 

-0.1 

9: 

0061 

1.5 

-0.1 

10: 

0061 

3.0 

-0.1 

11: 

0061 

4.5 

0.0 

12: 

0061 

0.0 

0.0 


FULL-PLANE CALCULATION 

AIRCRAFT SEMI -SPAN/DISK RADIUS: 0.60000E+01 M 

FLIGHT SPEED: 0.30000E+02 M/SEC 

NOMINAL RELEASE HEIGHT: 0.30000E+01 M 

VARIABLES FROM WAKE PLOT FILE: V W QQ 
SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

WAKE PLOT FILE TITLE: 

TEST WITH VELOCITIES AND QQ 
WAKE PLOT FILE MESH SIZES: (3,3) 

Y (M) : 0.0 0 . 2000E+02 0.4000E+02 

Z (M): 0.0 0 . 2000E+02 0.4000E+02 

WAKE PLOT FILE TIMES: 

0.0 SEC 

0 . 2000E+01 SEC 
0 . 5000E+01 SEC 
-0 . 1000E+01 SEC 
TOTAL NUMBER OF PARTICLES: 7 

DIAMETER: 0.10000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

$$$ WARNING: WAKE PLOT FILE INVOKED 
INITIAL TIME STEP: 0.15600E-01 SEC 

MAXIMUM TIME: 0.10000E+02 SEC 


TIME: 

0 . 6851E+00 

SEC 

TIME: 

0 . 1398E+01 

SEC 

TIME: 

0 . 2110E+01 

SEC 

TIME: 

0 . 2828E+01 

SEC 


INTEGRATION COMPLETE 

NASA AGDISP (MOD 2.0) PROGRAM END 


Figure 16a. AGDISP test case #5 terminal output 


NASA AGDISP (MOD 2.0) PROGRAM RESULTS 
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NASA AGDISP (MOD 2.0) PROGRAM RESULTS 


1: 

0000 

AGDISP TEST CASE # 

6 

2: 

0010 

10.0 1 


3: 

0020 

3 0 8.0 6.0 30. 

0 0 

4: 

0030 

15000.0 0.05 


5: 

0050 

0 1.0 30.0 


6: 

0060 

6 0 -1.0 100.0 

1.0 


HALF-PLANE CALCULATION 

AIRCRAFT SEMI -SPAN/DISK RADIUS: 0.80000E+01 M 

FLIGHT SPEED: 0.30000E+02 M/SEC 

NOMINAL RELEASE HEIGHT: 0.60000E+01 M 

HELICOPTER FORWARD ADVANCE: 0.50000E-01 

DOWNWASH VELOCITY: 0.53750E+01 M/SEC 

EFFECTIVE GAMMA: 0.12738E+01 M**2/SEC 

TURBULENCE FIXED VALUE: 0.10000E+01 (M/SEC)**2 

SCALE LENGTH MAXIMUM VALUE: 0.30000E+02 M 

TOTAL NUMBER OF PARTICLES: 6 

DIAMETER: 0.10000E+03 MICRONS 

SPECIFIC GRAVITY: 0.10000E+01 

INITIAL TIME STEP: 0.15600E-01 SEC 

MAXIMUM TIME: 0.10000E+02 SEC 

TIME: 0 . 6901E+00 SEC 

TIME: 0 . 1412E+01 SEC 

TIME: 0 . 2140E+01 SEC 

TIME: 0 . 2872E+01 SEC 

INTEGRATION COMPLETE 

NASA AGDISP (MOD 2.0) PROGRAM END 


Figure 17a. AGDISP test case #6 terminal output 
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APPENDICES 


The appendices contain variable definitions for the FORTRAN 
source code for AGDISP , AGPLOT and AGLINE. Explanation of the role 
of the subroutines may be found in Section 8 of this manual. 

The variables summarized in Appendices A, B and C do not, as 
a rule, include the following: 


1 . DO loop indices ; 

2. dummy variables in subroutine statements; 

3. temporary variables (holding the results of intermediate 
calculations); or 

4. variables in the scientific subroutines AGBZQ, AGBZT, 
AGLQD , AGLQS and AGRTF. 

Abbreviations have been used whenever possible; thus, deriv = 
derivative; eqn = equation; relax = relaxation; se = super- 
equilibrium; wrt = with respect to; etc. Vector variables are 
identified as such. Duplicated variable names occurring in AGPLOT 
are not repeated in AGLINE unless their function has changed. 



APPENDIX A: 


AGDISP Variable Summary 



NAME 


SUBR 


DESCRIPTION 


AGV AGBZG vector: absolute value of Betz circulation distribution 

APRP AGDISP initial frontal area of propeller 

AS AGDISP airplane planform area 

ASV AGP AC vector: integrand of average plant area density 

AV AGDISP vector: inputted plant area density, function of ZV 

AGPAC 
AGPAD 

AGWKI vector: wake plot file contents at a specified time 

AGWKR 

AGVKS 

B AGCOR ratio of particle relax time to turbulent relax time 

C AGCOR ratio of actual time to turbulent relax time 

CD AGDISP airplane drag coefficient 

AGPAC effective plant canopy vortex drag coefficient 

CPQ AGDISP constant for propeller turbulence evaluation 

AGPRP 

CPR AGDISP constant for propeller radius evaluation 

AGPRP 

CTA AGDISP cosine of terrain angle 

DCUT AGDEC particle diameter below which AGDISP terminates 

AGDISP 

DECAY AGDIF particle decay constant 

DELT AGBZI integration time step 

AGPAC 
AGVCH 

DELTA AGTUR spatial distance across which se derivs are computed 

DENF AGDEC particle specific gravity 

AGDISP 

DG AGBZG Betz fully rolled up vortex circulation strength 

AGBZI 

DGMM AGBZG error tolerance on computed Betz circulation derivs 

DGV AGBZD vector: Betz circulation derivs, functions of YV 

AGBZG 
AGBZI 
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NAME 


SUBR DESCRIPTION 


DIAM AGDEC particle diameter 

AGDISP 

DIST AGB2G initial vertical distance from ground to Betz vortex sheet 

AGDISP initial vertical distance from ground to tip vortex center 

DMAX AGSUP se maximum local spatial velocity deriv 

DMCV AGDISP vector: particle volume ratio at surface impact 

AGINT 

DNV AGBZI vector: new time derivs of Betz vortex eqns 

AGINT vector: new time derivs of particle eqns 

DOV AGBZG vector: old time derivs of Betz vortex eqns 

AGBZI 

AGINT vector: old time derivs of particle eqns 

DSYM AGBZG vector: unrolled Betz sheet length left-of -centroid 

AGBZI 
AGVEL 

DSYP AGBZG vector: unrolled Betz sheet length right-of -centroid 

AGBZI 
AGVEL 

DT AGBZI integration time step 

AGDISP 
AGINT 

DTAU AGCOR particle relax time 

AGDEC 
AGDIF 
AGDISP 

DTEMP AGDISP evaporation wet bulb temperature difference 

DTMN AGDIF minimum particle relax time, function of DIAM 

AGINT 

DTV AGBZG vector: Betz integration time step for each vortex 

AGBZI 

DU AGDEC particle velocity relative to local mean wind velocity 

DV AGBZD vector: new time derivs of Betz vortex eqns 

AGDIF vector: new time derivs of particle eqns 

AGMAT vector: se turbulence solution 
AGSUP 

DVDY AGSUP deriv of horizontal velocity wrt horizontal distance 

AGTUR 
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NAME 


SUBR DESCRIPTION 


DVDZ 

AGSUP 

AGTUR 

DWDY 

AGSUP 

AGTUR 

DWDZ 

AGSUP 

AGTUR 

DY 

AGBZG 

AGBZI 

DZ 

AGBZG 

AGBZI 

AGDISP 

DZBP 

AGDISP 

D1 

AGBZD 

D2 

AGBZD 

EPS 

AGDEC 

AGSUP 

AGTUR 

ETA 

AGDISP 

ETAU 

AGDEC 

AGDISP 

AGINT 

FACE 

AGBZG 

AGDISP 

AGPAC 

AGVEL 

FACR 

AGBZG 

AGDISP 

AGPAC 

AGVEL 

FL 

AGPAC 

AGPAD 

FR 

AGPAC 

AGPAD 

G 

AGSVE 

GAMMA 

AGDISP 


deriv of horizontal velocity wrt vertical distance 

deriv of vertical velocity wrt horizontal distance 

deriv of vertical velocity wrt vertical distance 

Betz fully rolled up vortex radius 
horizontal increment in Betz vortex location 

distance above DIST to biplane vertical location 
vertical increment in Betz vortex location 
incremental vertical distance 

incremental vertical distance from DIST to biplane 
Betz left-of-centroid spatial circulation deriv 
Betz right-of-centroid spatial circulation deriv 
particle decay constant 

se incremental turbulence level for root finding 
se continuity equation error 

airplane propeller efficiency 

particle evaporation relax time 

vector: plant canopy circulation reduction -Y vortices 

vector: plant canopy circulation reduction +Y vortices 

vector: plant area density integral -Y vortices 

vector: plant area density integral +Y vortices 

Betz unrolled vortex circulation strength 
triangular/rectangular vortex circulation strength 
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NAME 

SUBR 

DESCRIPTION 


GSAV 

AGBZG 

AGBZI 

AGDISP 

AGVEL 

vector: Betz average sheet circulation per unit 

length 

GV 

AGBZG 

AGBZI 

AGDISP 

vector: inputted Betz circulation, function of 

YV 

G1 

AGBZI 

Betz left-of-centroid circulation value 


G2 

AGBZI 

Betz right-of-centroid circulation value 


G2PI 

AGBZG 

vector: circulation divided by two pi (6.2831853) 


AGBZI 

AGDISP 

AGPAC 

AGVEL 



HHEL 

AGDISP 

AGINT 

AGVEL 

initial height of helicopter rotor plane above 

surface 

HTPAD 

AGDISP 

AGINT 

AGPAC 

AGTUR 

AGVEL 

plant canopy maximum height above surface 


I CARD 

AGDISP 

running count of input cards 


ICMX 

AGDISP 

maximum number of input cards 


ICV 

AGDISP 

vector : input card numbers 


IER 

AGMAT 

AGSUP 

se error flag 


IFL 

AGPAC 

AGPAD 

vector: -Y vortex within plant canopy flag 


IFR 

AGPAC 

AGPAD 

vector: +Y vortex within plant canopy flag 


IMAX 

AGSUP 

se maximum iteration limit 


I OUT 

AGINT 

AGSAV 

terminal output flag 


ISW 

AGINT 

vector: individual particle-to-ground impact flag 


AGSAV 




NAME 

SUBR 

DESCRIPTION 

IStfC 

AGINT 

AGSAV 

summed particle-to-ground impact flag 

IY 

AGWKI 

horizontal vector location within wake array data 

IZ 

AGWKI 

vertical vector location within wake array data 

LBP 

AGDISP 

biplane flag 

LCV 

AGDISP 

vector: admissible input card entries 

LENDF 

AGDEC 

AGDISP 

AGINT 

end-of -run flag 

LERF 

AGBZG 

Betz error flag 

LEVAP 

AGDEC 

AGDISP 

AGINT 

evaporation flag 

LHFPL 

AGDISP 

half-plane/full-plane flag 

LMCRS 

AGDISP 

AGVEL 

mean cross wind flag 

LMVEL 

AGDISP 

AGINT 

AGVCH 

AGVEL 

mean velocity flag 

LOC 

AGDISP 

input card counter 

LOCA 

AGDISP 

plant area density input card counter 

LOCB 

AGDISP 

Betz circulation input card counter 

LPART 

AGDISP 

particle flag 

LPRP 

AGDISP 

AGINT 

AGTUR 

AGVCH 

AGVEL 

propeller flag 

LQQSE 

AGDISP 

AGTUR 

turbulence flag 

LV 

AGDISP 

AGINT 

vector: input card flags 

vector: positive definite variables flag 

LZERO 

AGDISP 

centerline particle flag 
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NAME 

SUBR 

DESCRIPTION 

L10 

AGDISP 

card 0010 flag 

L20 

AGDXSP 

card 0020 flag 

L21 

AGDISP 

card 0021 flag 

L22 

AGDISP 

card 0022 flag 

L25 

AGDISP 

card 0025 flag 

L28 

AGDISP 

card 0028 flag 

L30 

AGDISP 

card 0030 flag 

L50 

AGDISP 

card 0050 flag 

L60 

AGDISP 

card 0060 flag 

L61 

AGDISP 

card 0061 flag 

L65 

AGDISP 

card 0065 flag 

ME 

AGBZD 

AGBZG 

AGBZI 

ending pointer for Betz vortex circulation 

MEV 

AGBZG 

AGBZI 

vector: ME values for all Betz vortices 

MS 

AGBZD 

AGBZG 

AGBZI 

starting pointer for Betz vortex circulation 

MSV 

AGBZG 

AGBZI 

vector: MS values for all Betz vortices 

MT 

AGBZD 

AGBZG 

AGBZI 

total discrete points in Betz vortex circulation 

MX 

AGBZD 

AGBZG 

AGBZI 

position of maximum Betz circulation deriv 

MXV 

AGBZG 

AGBZI 

vector: MX values for all Betz vortices 

NBTZ 

AGBZG 

AGBZI 

number of Betz vortices 

NDAT 

AGDISP 

card input file unit number 
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NAME 


SUBR DESCRIPTION 

NENDF AGWKR wake plot file record position counter 

AGWKS 

NEXTF AGWKI wake plot file extrapolation flag 

AGWKS 

NF AGWKR wake plot file record read counter 

NGAM AGBZG number of discrete inputted Betz circulation data points 

AGDISP 

NMAX AGWKS maximum admissible grid positions in wake plot file 

NOUT AGBZG terminal output unit number 

AGDISP 
AGMAT 
AGPAD 
AGSAV 
AGSUP 
AGWKS 

NPAD AGDISP number of discrete inputted plant area density data points 

AGPAC 
AGPAD 

NPLT AGDISP plot output file unit number 

AGSAV 

NPRT AGBZG printer output file unit number 

AGBZI 
AGDISP 
AGPAC 
AGSAV 
AGWKI 
AGWKR 

NSAV AGDISP frequency of writes to plot output file 

AGINT 

NV AGWKS vector: number of admissible variables in wake plot file 

NVAR AGDIF total number of particles 

AGDISP 
AGINT 
AGSAV 
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NAME 

SUBR 

NVOR 

AGBZG 

AGBZI 

AGDISP 

AGPAC 

AGTUR 

AGVCH 

AGVEL 

NWPF 

AGVKR 

AGWKS 

NWPV 

AGWKR 

AGWKS 

NY 

AGWKI 

AGWKR 

AGWKS 

NZ 

AGWKI 

AGWKR 

AGWKS 

PGBP 

AGDISP 

AGBZG 

AGBZI 

PSBP 

AGDISP 

AGBZG 

AGBZI 

QQ 

AGDIF 

AGTUR 

QQMX 

AGDISP 

AGTUR 

QQPRP 

AGDISP 

AGPRP 

AGTUR 

R 

AGTUR 

AGVEL 

REYNO 

AGDEC 

RHEL 

AGDISP 

AGVEL 

RLIM 

AGBZG 

AGDISP 

AGVEL 


DESCRIPTION 

total number of vortices 


wake plot file unit number 

total number of wake plot file variables 

total number of horizontal points in wake plot fil 

total number of vertical points in wake plot file 

biplane circulation factor 

biplane semispan factor 

local background turbulence level 
constant background turbulence level 
propeller turbulence level 

radius 

Reynolds number based on relative velocity 
helicopter blade radius 

limiting cutoff radius for triangular wing loading 



NAME 

SUBR 

DESCRIPTION 


RPRP 

AGDISP 

AGPRP 

AGTUR 

AGVEL 

radius of propeller flow field effect. 

function of 

S 

AGDISP 

AGSVE 

AGVEL 

airplane semispan 

length of unrolled Betz vortex sheet 


SL 

AGDIF 

AGTUR 

turbulent macroscale length 


SLMX 

AGDISP 

AGTUR 

maximum turbulent macroscale length 


SRV 

AGBZG 

AGBZI 

AGDISP 

AGPAC 

vector: vortex core radius 


STA 

AGDISP 

sine of terrain angle 


T 

AGCOR 

AGDEC 

AGDIF 

AGINT 

AGSAV 

AGWKR 

integration time 


TA 

AGDISP 

terrain slope angle 


TDOT 

AGDISP 

propeller shaft rpm 


TFE 

AGWKR 

wake plot file time just beyond desired 

time 

TFS 

AGWKR 

wake plot file time just before desired 

time 

TIME 

AGBZI 

AGDISP 

integration time 


TMAX 

AGDEC 

AGDISP 

AGINT 

maximum integration time 


TMCV 

AGDISP 

AGINT 

vector: time of particle-surface impact 


TR 

AGWKR 

AGWKS 

wake plot file interpolation time 


U 

AGDISP 

cross wind reference velocity 
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NAME 

SUBR 

DESCRIPTION 

UO 

AGDISP 

AGINT 

aircraft flight speed 

USK 

AGDISP 

AGVEL 

constant for cross wind evaluation 

UU 

AGSUP 

AGTUR 

local background axial turbulence component 

uv 

AGCOR 

AGDIF 

particle velocity/background velocity correlation 

ux 

AGCOR 

AGDIF 

particle position/background velocity correlation 

UY 

AGMAT 

AGSUP 

deriv of axial velocity wrt horizontal distance 

UZ 

AGMAT 

AGSUP 

deriv of axial velocity wrt vertical distance 

V 

AGDIF 

AGDISP 

AGSVE 

AGVEL 

horizontal velocity 

VPRP 

AGDISP 

AGPRP 

AGVEL 

propeller tip velocity 

W 

AGDIF 

AGSUP 

AGTUR 

local background horizontal turbulence component 

VY 

AGMAT 

AGSUP 

deriv of horizontal velocity wrt horizontal distance 

VZ 

AGMAT 

AGSUP 

deriv of horizontal velocity wrt vertical distance 

V 

AGDIF 

AGDISP 

AGSVE 

AGVEL 

vertical velocity 

WHEL 

AGDISP 

AGINT 

AGVEL 

helicopter downwash velocity at blade plane 

WT 

AGDISP 

helicopter weight 
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NAME 

SUBR 

DESCRIPTION 

WTAU 

AGCOR 

AGDIF 

turbulent relax time 

WW 

AGDIF 

AGSUP 

AGTUR 

local background vertical turbulence component 

WY 

AGMAT 

AGSUP 

deriv of vertical velocity wrt horizontal distance 

WZ 

AGMAT 

AGSUP 

deriv of vertical velocity wrt vertical distance 

X 

AGPRP 

AGSVE 

AGWKI 

axial distance measured from particle release point 

XMU 

AGDISP 

helicopter forward advance ratio 

XNV 

AGBZI 

AGINT 

vector: new time values of Betz vortex eqns 
vector: new time values of particle eqns 

XO 

AGINT 

AGVEL 

axial distance measured from particle release point 

XOV 

AGBZG 

AGBZI 

AGDISP 

AGINT 

vector: old time values of Betz vortex eqns 
vector: old time values of particle eqns 

XPR 

AGDISP 

AGPRP 

propeller axial virtual origin 

XV 

AGBZD 

AGDIF 

AGSAV 

vector: new time values of Betz vortex eqns 
vector: new time values of particle eqns 

XI 

AGBZI 

Betz integral for left-of -centroid 

X2 

AGBZI 

Betz integral for right-of -centroid 

Y 

AGSVE 

AGTUR 

AGVEL 

AGWKI 

horizontal distance 

YBAL 

AGBZG 

AGBZI 

AGDISP 

AGTUR 

AGVCH 

AGVEL 

vector: horizontal location of 2nd quadrant vortices 
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NAME 

SUBR 

YBAR 

AGBZG 

AGBZI 

AGDISP 

AGTUR 

AGVCH 

AGVEL 

YE 

AGBZI 

YHEL 

AGDISP 

AGVCH 

AGVEL 

YNV 

AGBZI 

YOV 

AGBZG 

AGBZI 

YPRP 

AGDISP 

AGTUR 

AGVCH 

AGVEL 

YS 

AGBZI 

YV 

AGBZD 

AGBZG 

AGBZI 

AGDISP 

AGWKI 

AGWKS 

Y1 

AGBZD 

AGBZI 

Y2 

AGBZD 

AGBZI 

Z 

AGDISP 

AGSVE 

AGTUR 

AGVEL 

AGWKI 

ZBAL 

AGBZG 

AGBZI 

AGDISP 

AGPAC 

AGTUR 

AGVCH 

AGVEL 


DESCRIPTION 

vector: horizontal location of 1st quadrant vortices 


right -of -centroid maximum location of Betz rollup 
horizontal location of helicopter centerline 

new horizontal location of Betz vortex centerline 
vector: old horizontal location of Betz vortex centerlines 

horizontal location of propeller centerline 


left-of-centroid minimum location of Betz rollup 
vector: inputted discrete wing stations for Betz GV 

vector: horizontal spatial locations in wake plot file 
left-of-centroid location of Betz rollup point 
right-of -centroid location of Betz rollup point 
vertical distance 


vector: vertical location of 2nd quadrant vortices 
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NAME 


SUBR DESCRIPTION 


ZBAR AGBZG vector: vertical location of 1st quadrant vortices 

AGBZI 
AGDISP 
AGPAC 
AGTUR 
AGVCH 
AGVEL 

ZHEL AGDISP vertical location of helicopter centerline 

AGINT 
AGVCH 
AGVEL 

ZNV AGBZI new vertical location of Betz vortex centerline 

ZO AGDISP surface roughness height 

AGVEL 

ZOPAD AGDISP plant canopy effective surface roughness height 

AGVEL 

ZOV AGBZG vector: old vertical location of Betz vortex centerlines 

AGBZI 

ZPRP AGDISP vertical location of propeller centerline 

AGTUR 
AGVCH 
AGVEL 

ZV AGDISP vector: inputted vertical locations of plant density AV 

AGPAC 
AGPAD 

AGWKI vector: vertical spatial locations in wake plot file 
AGWKS 
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APPENDIX B: 


AGPLOT Variable Summary 



NAME SUBR DESCRIPTION 

ANSWR AGPLOT interactive user response 

AGPLT 

AV AGPLT vector: particle solution from plot output file 

AGS PD 

CT AGSPD particle direction cosine 

DELT AGPLOT tag time increment 

AGPLT 

DIST AGPLOT initial vertical distance from surface to aircraft wing 

AGSPD 

DMCV AGPLOT vector: particle volume ratio at surface impact 

DV AGGRD vector: particle volume ratio at surface impact 

AGPLOT vector: horizontal spread and location at surface impact 

DZBP AGPLOT incremental vertical distance from DIST to biplane 

GDLT AGGRD ground deposition plot scale interval 

GMAX AGGRD ground deposition plot scale maximum 

GMIN AGGRD ground deposition plot scale minimum 

GV AGGRD vector: ground deposition 

ICARD AGPLOT running count of AGDISP input cards 

ICV AGPLOT vector: particle-surface impact flag 

IE AGPLT number of scale tick marks 

IV AGGRD vector: particle-surface impact flag 

IX AGCHK number of scale divisions 

IXD AGDRW horizontal screen interval size equivalent 

AGPLT 

IXE AGDRW horizontal screen maximum size equivalent 

AGPLT 

IXN AGDRW vector: horizontal screen coordinate equivalent 

IXS AGDRW horizontal screen minimum size equivalent 

AGPLT 

IXT AGPLT horizontal screen scale divisions 
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NAME 


SUBR 


DESCRIPTION 


IXV 

AGPLT 

IXZ 

AGPLT 

1X1 

AGPLT 

1X2 

AGPLT 

1X3 

AGPLT 

I YD 

AGDRW 

AGPLT 

IYE 

AGDRW 

AGPLT 

IYN 

AGDRW 

IYS 

AGPLT 

IYT 

AGPLT 

IYV 

AGPLT 

IYZ 

AGPLT 

IY1 

AGPLT 

IY2 

AGPLT 

IY3 

AGPLT 

LB OX 

AGPLT 

LENDF 

AGPLT 

LFOLD 

AGGRD 

AGPLOT 

AGPLT 

LHF 

AGGRD 

LHFPL 

AGPLOT 

LMPLT 

AGPLOT 

AGPLT 

LSC 

AGPLT 

LSPLT 

AGPLOT 

AGPLT 


vector: terminal horizontal screen limits 
horizontal screen zero axis coordinate equivalent 
horizontal screen pointer for mean particle index 
horizontal screen pointer for right-of -particle index 
horizontal screen pointer for left-of-particle index 
vertical screen interval size equivalent 

vertical screen maximum size equivalent 

vector: vertical screen coordinate equivalent 

vertical screen minimum size equivalent 

vertical screen scale divisions 

vector: terminal vertical screen limits 

vertical screen zero axis coordinate equivalent 

vertical screen pointer for mean particle index 

vertical screen pointer for right-of -particle index 

vertical screen pointer for left-of-particle index 

tag time pointer 

end-of-plot flag 

symmetry plot flag 

half -plane/full-plane flag 
half-plane/full-plane flag 
mean particle trajectory plot flag 

equal plot scales flag 

standard deviation trajectory plot flag 
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NAME 

SUBR 

DESCRIPTION 

LTERM 

AGGRD 

AGMDL 

AGPLOT 

AGPLT 

terminal type flag 

LTPLT 

AGPLOT 

AGPLT 

tag plot flag 

LZERO 

AGPLOT 

centerline particle flag 

LZR 

AGGRD 

centerline particle flag 

NINU 

AGGRD 

AGPLOT 

AGPLT 

terminal input unit number 

NOUT 

AGCHK 

AGGRD 

AGMDL 

AGPLOT 

AGPLT 

terminal output unit number 

NPLT 

AGPLOT 

AGPLT 

plot input file unit number 

NPTS 

AGDRW 

AGGRD 

total number of plotted particles 

NTOT 

AGGRD 

ground deposition discretization 

NVAR 

AGGRD 

AGPLOT 

AGPLT 

total number of particles 

PSBP 

AGPLOT 

biplane semispan factor 

RSPD 

AGSPD 

average standard deviation 

S 

AGPLOT 

aircraft semispan 

SRTPI 

AGGRD 

constant for ground deposition normalization 

SSS 

AGDRW 

AGGRD 

AGPLT 

symmetry parameter 

ST 

AGSPD 

particle direction sine 

SY 

AGGRD 

horizontal particle turbulent spread 

TA 

AGPLOT 

AGPLT 

terrain slope angle 
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NAME 

SUBR 

DESCRIPTION 


TIME 

AGPLOT 

AGPLT 

simulation time 


TMAX 

AGPLOT 

maximum simulation time 


TV 

AGPLT 

vector: simulation time 


XD 

AGDRW 

AGPLT 

horizontal screen interval value 


XDC 

AGPLT 

horizontal screen interval value for equal 

axes 

XDD 

AGPLT 

horizontal screen interval value 


XDLT 

AGCHK 

plot scale interval 


XMAX 

AGCHK 

plot scale maximum 


XMIN 

AGCHK 

plot scale minimum 


XS 

AGDRW 

AGPLT 

horizontal screen starting value 


XSS 

AGPLT 

horizontal screen starting value 


XTT 

AGPLT 

horizontal screen total value 


XV 

AGDRW 

AGGRD 

vector: horizontal screen coordinate 

vector: horizontal spread and location at surface impact 

Y 

AGGRD 

horizontal distance 


YD 

AGDRW 

AGPLT 

vertical screen interval value 


YDC 

AGPLT 

vertical screen interval value for equal axes 

YDD 

AGGRD 

AGPLT 

vertical screen interval value 


YDLT 

AGPLOT 

horizontal axis plot scale interval 


YE 

AGGRD 

ending horizontal distance along surface 


YM 

AGPLOT 

particle horizontal mean position from plot 

file 

YMAX 

AGPLOT 

horizontal axis plot scale maximum 


YMIN 

AGPLOT 

horizontal axis plot scale minimum 


YMMN 

AGPLOT 

minimum horizontal value for mean particle 

trajectories 
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NAME 

SUBR 

YMMX 

AGPLOT 

YP 

AGPLOT 

YS 

AGDRW 

AGGRD 

AGPLT 

YSMN 

AGPLOT 

YSMX 

AGPLOT 

YSS 

AGGRD 

AGPLT 

YTT 

AGGRD 

AGPLT 

YV 

AGDRW 

AGGRD 

AGPLT 

ZDLT 

AGPLOT 

ZM 

AGPLOT 

ZMAX 

AGPLOT 

ZMIN 

AGPLOT 

ZMMN 

AGPLOT 

ZMMX 

AGPLOT 

ZP 

AGPLOT 

ZSMN 

AGPLOT 

ZSMX 

AGPLOT 

zv 

AGPLT 


DESCRIPTION 

maximum horizontal value for mean particle trajectories 

particle horizontal standard deviation from plot file 

vertical screen starting value 

starting horizontal distance along surface 

vertical screen starting value 

minimum horizontal value for particle spread trajectories 
maximum horizontal value for particle spread trajectories 
vertical screen starting value 

vertical screen total value 

vector: vertical screen coordinate 
vector: horizontal distance 

vertical axis plot scale interval 
particle vertical mean position from plot file 
vertical axis plot scale maximum 
vertical axis plot scale minimum 

minimum vertical value for mean particle trajectories 
maximum vertical value for mean particle trajectories 
particle vertical standard deviation from plot file 
minimum vertical value for particle spread trajectories 
maximum vertical value for particle spread trajectories 
vector: vertical distance 
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APPENDIX C: 


AGLINE Variable Summary 


NAME 

SUBR 

DESCRIPTION 

CBAR 

AGLINE 

equivalent Gaussian magnitude 

CV 

AGELP 

vector: isopleth position value 

DFOM 

AGLINE 

desired f igure-of-merit 

DIAM 

AGEQD 

particle diameter 

ETAU 

AGEQD 

evaporation time constant 

FOM 

AGLINE 

equivalent Gaussian f igure-of-merit 

IGRDV 

AGLINE 

vector: ground interaction pointer 

I ZERO 

AGLINE 

number of zero spread points 

NPTS 

AGLINE 

desired sequence number 

SYB 

AGLINE 

equivalent Gaussian horizontal standard deviation 

SYV 

AGLINE 

vector: horizontal standard deviation 

SZB 

AGLINE 

equivalent Gaussian vertical standard deviation 

SZV 

AGLINE 

vector: vertical standard deviation 

XC 

AGEQD 

sum of concentration 

XFn 

AGEQD 

f igure-of-merit integrals (n=l,2,3,4) 

XY 

AGEQD 

sum of horizontal position 

XYS 

AGEQD 

sum of horizontal variance 

XZ 

AGEQD 

sum of vertical position 

XZS 

AGEQD 

sum of vertical variance 

YB 

AGLINE 

equivalent Gaussian horizontal position 

YV 

AGLINE 

vector: horizontal particle position 

ZB 

AGLINE 

equivalent Gaussian vertical position 
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